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Abstract 



We give an introduction to the recently established connection between su- 
persymmetric gauge theories and matrix models. We begin by reviewing 
previous material that is required in order to follow the latest developments. 
This includes the superfield formulation of gauge theories, holomorphy, the 
chiral ring, the Konishi anomaly and the large ./V limit. We then present 
both the diagrammatic proof of the connection and the one based on the 
anomaly. Our discussion is entirely field theoretical and self contained. 
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1 Introduction 

The scope of this review is to give a pedagogical introduction to some new 
techniques developed in the context of supersymmetric gauge theories, en- 
abling one to compute systematically non-perturbative quantities in the low 
energy effective physics. These new developments were initiated by a con- 
jecture put forward by Dijkgraaf and Vafa P where it was argued that 
non-perturbative quantities such as the glueball superpotential can be com- 
puted by means of a matrix model. 

The main motivation for the study of supersymmetric (SUSY) gauge 
theories is of course that they are expected to be directly relevant to the 
next round of high energy experiments, where evidence for the supersym- 
metric partners of the standard model particles should hopefully be found. 
To be phenomenologically viable, supersymmetry must of course be bro- 
ken in some way. While this is the main phenomenological context for the 
study of supersymmetry, there are further theoretical reasons for studying 
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SUSY gauge theories. These theories, even when preserving supersymme- 
try (to be precise, minimal J\f = 1 SUSY), are believed to share many key 
dynamical features with QCD such as confinement, dynamical generation 
of a mass gap in the gauge sector, chiral symmetry breaking and the re- 
lated non-trivial vacuum structure. All these features are non-perturbative 
in nature and notoriously difficult to address in real life QCD. With unbro- 
ken SUSY there are a number of powerful tools at one's disposal, such as 
non-renormalization theorems and holomorphy. While important dynamical 
issues such as confinement remain out of reach and have to be assumed, it is 
possible to make exact statements concerning the vacuum structure of the 
theory and the value of the vacuum condensates. 

The quantitative investigation of the non-perturbative aspects of SUSY 
gauge theories began in the early 80's, as reviewed in [2], and received new 
impetus during the 90's, see P]. Subsequently, in the second part of the 
90's it was realized that there are a huge number of possibilities to embed 
gauge theories within string theory, using D-branes and dualities. This 
led to an enormous amount of work in many directions, one of which 
culminated in the proposal of [Q. This conjecture, originally motivated 
in the above string theory context, was subsequently given a purely field 
theoretical interpretation in various different approaches. In |5] the emphasis 
was on the relation with the Seiberg-Witten solution of M = 2 pure gauge 
theory jHj, broken to N = 1 by the presence of a tree level superpotential 
for the adjoint matter field. In a purely perturbative argument was 
given, in which the glueball superfield was treated as a classical background 
field. Finally, in jSj the connection with a matrix model was made using 
generalized Konishi anomaly relations 

In this review we will concentrate on the purely field theoretical J\f=l 
approach to these new techniques, very much in the spirit of [7j and |S]. This 
choice is partly dictated by the wish to provide a self contained introduc- 
tion using the least amount of advanced results. Hopefully this will equip 
the reader with the basic knowledge to address the more advanced current 
literature on the subject. In this spirit, we have made some effort to discuss 
earlier results in J\f = 1 gauge theories to put the new developments into 
context. 

The review is organized as follows. In Section 2 we briefly present the su- 
perspace formulation of N = 1 gauge theory. This is the language in which 
all the results will be derived. We introduce chiral and vector superfields 
and show how they combine into the most general gauge theory with mat- 
ter. Section 3 discusses some basic dynamical properties of these theories 
such as the perturbative non-renormalization theorem, the holomorphy of 
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the Wilsonian superpotential, the linearity principle and the formulation of 
the low energy physics in terms of the glueball superfield. The review of the 
latest developments begins in Section 4 where the conjecture of Dijkgraaf 
and Vafa is presented and derived in a perturbative context. This sec- 
tion describes the computation of 7 by considering the simplest example, 
namely that of a cubic superpotential for matter in the adjoint representa- 
tion of U(N). We show explicitly the equivalence between the field theory 
computation and the matrix model concentrating on a simple diagram. Fur- 
thermore, we show that only planar diagrams are relevant to the problem. 
Section 5 introduces the concept of the chiral ring, which provides a useful 
tool to analyze the vacuum structure of the theory. To do so, we define chiral 
operators and derive their basic properties. We then introduce an equiva- 
lence relation between chiral operators leading to the concept of the chiral 
ring and discuss its quantum corrections. Section 6 describes the Konishi 
anomaly |Hj and its generalizations |Hj. Finally, in Section 7 we show how 
one can determine the effective superpotential, using the results of the two 
previous sections, and how this is related to the matrix model. This last 
section is closely based on (HJ. 

Appendix A fixes the notation. Appendix B reviews the one cut solution 
to the matrix model of ; needed for the complete solution of the example 
in Section 4. In Appendix C we make a humble attempt at reviewing the 
literature that uses the new techniques discussed here. With this we hope to 
partly correct for not having touched many interesting new developments in 
the main text. A second apology is due for omitting many important results 
in SUSY gauge theories previous to the Dijkgraaf and Vafa conjecture, let 
alone for not mentioning extended supersymmetry. Fortunately there are 
many extensive reviews covering these subjects. 

2 Superspace formulation of gauge theories 

The purpose of this section is to familiarize the reader with the techniques of 
superspace which has the advantage of making SUSY manifest, see e.g. 
for a list of references to the original papers. We begin by showing what 
a SUSY Lagrangian looks like in a non manifestly SUSY notation (the so 
called component notation) and then explain how it can be rewritten in a 
way so as to make SUSY manifest. The reader familiar with the formalism 
can skip directly to Section 3. 

The two typical examples of four dimensional SUSY Lagrangians are 
the so called Wess-Zumino (WZ) model (no gauge interactions, only matter 
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fields) and the pure gauge theory (only gauge interaction, no matter fields). 
The WZ model can be written as |12j : 

C = d^d^cj) — i^a^d^ip + ff + interactions, (1) 

where <p is a complex scalar, ip a Weyl fermion 1 and / an auxiliary field (a 
field without a kinetic term) that can be integrated out by its own (algebraic) 
equation of motion. The interaction terms can also be written explicitly 
and we will discuss them at length in the following. For the purpose of 
introducing the SUSY transformations however the kinetic terms (^Q) are 
enough. 

The pure gauge theory on the other hand has a Lagrangian JH] 

C = tr (-If^F^ - 2f\^D^\ + D 2 ^j , (2) 

where F^ u stands for the usual field strength, A is a Weyl fermion (the 
"gluino" or "gaugino") and D is, again, an auxiliary field. All fields in Q 
are matrix valued in some Lie algebra and the trace is the trace over the 
generators: tr (T a T b ) = \5 ab . One could also add a so called "topological 
term" proportional to e^po-tr F^ v F pa or, for U{1) factors in the gauge group, 
a so called Fayet-Iliopoulos (FI) term proportional to D but otherwise gauge 
invariance prevents any other type of renormalizable interaction. 

The most general four dimensional gauge theory with M = 1 SUSY is, 
loosely speaking, based on a combination of (0) and © as we are going to 
show. 

The SUSY transformations can be written explicitly as transformations 
on the fields. In the context of the WZ model, one can show that 

5(p = £ip 

5V = iv^d^ + U (3) 
5f = i^S M V 

in terms of one 2 constant Grassmann (i.e. anticommuting) Weyl spinor £ a 
(£« = (Co)*)) changes ^ by at most a total derivative, thus leaving the 
action S = f d 4 x£ invariant. A similar transformation can be shown to 

1 Weyl spinors are ubiquitous in SUSY Lagrangians, and the conventions used are 
slightly different from the more familiar ones involving Dirac spinors. We review the 
notation in Appendix A. 

2 Thus the name Af = 1. 
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leave (J2J) invariant, again up to a total spacetime derivative: 



6A» = _L(_iA^ + i^A) 
v2 



1 
1 



a = —((7^^ + ^) ( 4 ) 



SD = -^{-^D^X-D^O- 

While the invariance of Q under Q is straightforward, in checking the 
invariance of (J2J) under Q a more involved four Fermi term arises from the 
variation of the gauge potential inside D^, which can be shown to vanish 
using Fierz identities. For later convenience, in (jSJ) and (@J) we rescaled the 
SUSY parameter by a factor of 1/V2 from that which is commonly used. 

The formulation briefly sketched above is manifestly Lorentz and gauge 
invariant, but lacks manifest supersymmetry invariance. If we had such a 
manifest SUSY formulation we would not need to go through the explicit 
computations discussed in the previous paragraph. This task is accom- 
plished by noticing that the SUSY transformations can be written in terms 
of "translations in superspace" |14j . To understand what that means, let us 
consider an ordinary translation of a scalar field 4>(x). We know that there 
exist an operator V 11 such that 

(j){x + a) = e- iaV <j){x)e iaV . (5) 

When a is small we can expand both sides. On the one hand: 

<j){x + a) = (l-m^ + . . .)(j){x)(l + ia^V^ + . . .) = <j>(x) -ia^[V^ <f>(x)] + . . . 

(6) 

on the other hand 

4>(x + a) = (j)(x) + a^d^(p{x) + . . . (7) 

Thus: 

[V^^(x)]=id^(x) = P^(x). (8) 

Hence, we have two objects playing the role of "momentum": which is 
the generator of translations and = id^ which gives the representation 
of Vp, in field space. These two objects are often written using the same 
symbol. 

Supersymmetry is the generalization of this to "superspace", that is 
space parameterized by the usual coordinates and new complex Grass- 
mann (anti-commuting) coordinates 9 a and 9 a = (9 a )* . A "superfield" is a 
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function in superspace: Y(x, 9, 9). Since (9\) 2 = (#2) 2 = and similarly for 
9 due to anticommutativity, the Taylor expansion of Y in 9 and 9 terminates: 

Y{x,9,9) = (j)(x)+9T](x) + 9x(x) + 9 2 m(x) + 9 2 n(x)+9a^9A^x) + 
9 2 9\{x) + 9 2 9^{x) + 9 2 9 2 d(x). (9) 

We have not indicated it explicitly but Y could carry some spin indices that 
would reflect in the component structure. The above formula refers to a 
Lorentz scalar for simplicity. 

In analogy with usual translations, we now define translations in super- 
space by a quantity (a, £, £) as: 

9 a ► 9 a + £a 

h - 0d + £d (10) 
-»• a^ + e^ + ^a^- ^ct^, 

where one should note the crucial addition of the fermionic bilinears to 
the shift of x^. They are of course needed so that the composition of two 
fermionic translations yields a translation in x M . 

The operators generating such transformation on the superfields are V, 
Q and Q obeying 

Y(x" + a»+ % -9o»t - l -^9, 9 a + 9 a + h) = 
e- iar+ t Q -&Y(x>*, 9 a , h)e iaV -^+~&. (11) 

Not surprisingly, the operators Q and Q are also represented by differential 
operators acting on superfields, as one can check: 

[Qa,Y] = Q a Y, [Q a ,Y]=QaY, (12) 

where on the left hand side we have a commutator or anticommutator de- 
pending on the spin of Y . Explicitly: 

Q a = d a - fa. 9%, Q a = 8 a - (13) 

One can show that P, Q and Q obey the SUSY algebra 

{Qa,Qd} = -<A- (14) 
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Having seen that SUSY is a translation in superspace, any integral over 
superspace of a superfield Y: 

J d 4 xd 2 ea 2 e y, (15) 

will be manifestly supersymmetric because the integration measure is trans- 
lational invariant by construction (see Appendix A). Of course, Y has to be 
a scalar superfield and can arise as the composition of other superfields. 

In the case of translations we were done at this point. Here we have 
to perform an extra step because the generic superfield Y above forms a 
reducible representation of the SUSY algebra (|14j) . that is, it is possible to 
restrict the form of Y to contain only a subset of the original fields that are 
still mapped into each other by SUSY. It should be obvious that Y is "too 
big" because we would like, for instance, to formulate the WZ model (|T|) 
which does not have any vector fields or the pure gauge theory (J2J which 
does not have any scalars. The two ways of reducing the field content in Y, 
that are going to be used are: (i) The chiral (or antichiral) projection and 
(ii) The real projection. 

2.1 Chiral superfields 

A chiral superfield can be constructed by noticing that there exist two dif- 
ferential operators: 

D Q = d a + Z -a»J%, Ai = B a + (16) 

that anticommute with all the supercharges (|13j) . This means that if Y is a 
superfield, so are D a Y and D a Y and it is consistent with SUSY to set one 
of them to zero. (Setting both would result in requiring Y = const.) We call 
chiral superfield a superfield $ obeying = 0. (Analogously, D a <l> = 
defines an antichiral superfield.) 

The solution to = can be obtained by observing that the (com- 
plex) combination y^ = x^ + ^Oa^O obeys D&y 11 = 0, and so does trivially 
9 a - Thus 

^ = ( j ) (y) + e^{y)+e 2 f(y). (17) 

is the general expression for a chiral superfield. Expanding (|17[) around x 
yields 

$ = (j){x) + l -Qa^Qd^{x) - ^6 2 6 2 Ucj)(x) + 6^(x) - l -e 2 d^{x)^6 + 6 2 f{x) 
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Let us begin by considering the action of Q on <I>. Recalling that SUSY 
can be written as a translation in superspace, we have, on the one hand 3 : 

^* = £Q* = ^fo)+£0/(y), (19) 

on the other hand: 

= [£Q, $] = [£Q, 0] + Q, ^] + 2 [£Q, /]. (20) 

Thus, the transformation properties © are reproduced, now written in 
terms of the SUSY generators 

[Qa,4>] = i'oc 

{Qa,M = ~e a pf (21) 
[QaJ] = 0. 

Similarly, one can check that 

[QaA] = o 
{Q&M = -™%cP^ (22) 

In order to write a SUSY action we recall that the superfield Y in (|15l) 
can be the composition of other superfields. For instance, given a chiral 
superfield $ one can write 

J d 4 xd 2 ^d 2 ^ $$, (23) 

which has been chosen as an example because by expanding and $ and 
integrating over 6 and 9 we obtain precisely the kinetic terms in the WZ 
model (UJ). 

It is often stated that there is an exception to (|15j) when the integrand 
happens to be a chiral superfield E (or a product of such). In this case, 
the integral over the whole superspace vanishes, however, the integral over 
chiral superspace 

J d 4 xd 2 # S (24) 

is still manifestly SUSY (translation invariant) as can readily be seen in the 
formalism where x^ 1 is replaced by so that 6 never appears. Equivalently, 

3 The SUSY transformation of the coordinate y M is independent of £. We take formally 
£ = since we are only interested in Q and not Q. 
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it can be noticed from eqs. (|21[) and (|22|) that the highest component / of 
a chiral superfield (which is singled out by the 9 2 integral) transforms as a 
total spacetime derivative. 

However, it should be appreciated that (|24|) is actually more general 
than (|15|) in the sense that any integral of type (|15|) can be written as 
()24j) . Remember that integrals over a Grassmann variable are equivalent to 
derivatives and that total spacetime derivatives do not contribute to the 
bosonic integral 4 so that one can write 



The superfield D 2 Y is manifestly chiral, regardless of the form of Y because 
three always give zero since they anticommute and there are only two 
possible values a can take. We call such terms "chirally exact". On the 
other hand, not all terms of the form (|24j) can be written as tjl5|) in terms 
of a local integral. For instance, the quantity 



for some chiral superfield <3?, cannot be expressed as (|15|) because we do not 
have any derivative that can be reconverted into a d 2 9. We call terms such 
as (|26|) that cannot be written as local integrals over the whole superspace 
"F-terms" and all the others, including (|25|) . "D-terms". 

The above expression (|2l-{jl is so compact that it may be a bit confusing at 
first - for instance, if one tried to obtain the equations of motion by varying 
<& one would get the nonsensical result $ = 0. The error, of course, is that 
the field $ obeys = but the action (f2*5|) doesn't know that. An easy 
way to enforce this is to use the above tricks and write 



Now the fields are unconstrained in (anti)-chiral superspace and varying 
with respect to <£> yields D 2 <& = which implies the free equations of motion 
for (f>, ip and /. 

Interactions can be introduced in the same way. Since the simplest non- 
vanishing integral over all of superspace ()23|) already involves derivatives, we 
must turn to terms like (|26[) . In fact any holomorphic function H^( ( I ) ) (i.e. 

4 We need not worry about boundary terms at this stage. 




(25) 




(26) 




(27) 
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such that dW/d$> = 0) to be called the superpotential, will be acceptable 
in the sense that if <I> is chiral, so is W(3>): 



, . 8W(<!>) - 8W($) - - 
D«W($) = —q^-D** + —g^D^ = 0. (28) 

Thus 

J d*xd 2 9 W($) + J d 4 xd 2 W($) (29) 

gives a super symmetric interaction term that can be added to the kinetic 
terms in (|23[1 . If renormalizability is an issue (and it will not be in many 
following cases since we will be dealing with effective theories) we must 
restrict W to be at most a cubic polynomial. 



2.2 Vector superfields 

To introduce gauge interactions we need a second type of superfield - the 
vector superfield. This is just the original Y (conventionally called V in this 
case and Lie algebra valued) on which we impose the reality condition 

V = V ] . (30) 

This condition is trivially supersymmetrically covariant because if V is the 
most arbitrary superfield, so is and equating the two is a superfield 
equation. In components, displaying the Lie algebra index: 

V a = C a {x) + i9 X a (x) - idx a {x) + i9 2 m a {x) - id 2 m a {x) - BaNA^x) 
+V2i0 2 9(\ a (x) + -^a^d^ix)) ~ V2i0 2 B(\ a {x) 
+ ^a^x a (x)) + 9 2 9 2 (D a (x) - \aC a (x)), (31) 

where C a , A^ and D a are real and the shifts in A a and D a by derivatives 
of x a an d C a are for later convenience. If T a are the generators of the Lie 
algebra ([T a ,T b ] = if abc T c ), one can write a matrix valued object V = V a T a 
obeying V = V' because V a are real and T a hermitian. 

Let us now see the relation between the vector superfield and gauge 
invariance. Consider matter described by a chiral superfield $ transforming 
in some representation (reducible or irreducible) of the Lie algebra. The 
usual gauge transformation 

$ _» $' = e ~ iA $ (32) 
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keeps this form but now A = A a T a must also be promoted to a chiral 
superfield because otherwise would no longer be chiral. However there is 
a problem with the kinetic term because, being A complex 

_> <j> e iA e ~ iA $ (33) 

Thus, we must compensate by inserting e v between $ and $ and postulate 
the transformation law: 

e V _ e V> = e -iA e V e iA (34) 

so that ^e^^ is gauge invariant. 

At first sight, it is not obvious that (|34|) implies an infinitesimal gauge 
transformation like 

SAfj, = d^a + ilA^a) (35) 

but, in fact, (|35[) is included in the transformation (|34|) of the components 
of V as one should check at least for an Abelian group, where (|34|) simplifies 
to V V = V - iA + iA. 

Using Q34j) it is possible to show that one can gauge away many of the 
fields in the original definition of V and go to the so called Wess-Zumino 
gauge where 

v a = -ea^eA^x) + V2i9 2 ex a { x ) - V2i9 2 ex a {x) + e 2 e 2 D a (x). (36) 

In this form, V is no longer a superfield - this is not a SUSY gauge. Indeed 
acting with Q or Q would regenerate the terms that have been set to zero in 
the WZ gauge. However, these new terms can be once again removed by a 
new (field dependent) gauge transformation. In other words, in this gauge 
the SUSY algebra closes up to a gauge transformation. 

As for the superpotential, no modification is required. As long as W is 
holomorphic and gauge invariant in the usual sense, the gauge parameter 
can be promoted to a full chiral superfield because its complex conjugate is 
never used. This is often stated by saying that the F-terms are invariant 
under the complexified gauge group. 

We now turn to finding the kinetic term for V, i.e. the manifestly SUSY 
way of writing ((21). From the transformation law (|34|) and the fact that 
appears explicitly in V it is clear that we must act with some derivative to 
construct an object to be identified with the supersymmetric field strength. 

It turns out that the right expression is 

W a = D 2 {e- y (D a e v )) (37) 
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where we have defined the manifestly chiral superfield W a 5 . The antichiral 
superfield Wa has a similar definition, obtained from W a = (W a )* . 

It is a nice exercise to check that, under a gauge transformation (|34|) . 
W a transforms as 

W a ^W' a = e~ iA W a e iA . (38) 

Note the difference between (|38|) and (j34j) - in (|38[) there is no A appearing 
and so it is possible to make a gauge invariant combination e.g. by taking 
traces. In checking ()38(l one has to notice that, first (|34|) implies e~ v — > 
e -iA e - ^e iA , second that the A terms commute through D a and cancel and, 
finally, that the spurious term D 2 D a (e lA ) vanishes. 

W a can be easily computed in the WZ gauge to give 

W a = -V2i\ a (y) + 6 a D{y) - ia^%F^(y) + ^/.^A^), (39) 

where, as usual, expanding in y would introduce additional spacetime deriva- 
tives. Since W a starts off with the gaugino, it is often called the gaugino (or 
gluino) superfield. 

The pure gauge theory action (jJJ) can be written in a manifestly SUSY 
way using the field W a : 

J d A xd 2 8trW a W a + J d 4 xd 2 6 tr W&W*. (40) 

For later purpose, it is convenient to reformulate expression (|37|) for the 
field strength W a in terms of (anti)commutators of some new "covariant" 
derivatives. Define 

Va = A* and V a = e~ v D a e v , (41) 

the second equation being intended as an operator, that is D a acts on e v 
and whatever follows it. Given the two new derivatives, one can define 
"covariantly chiral/antichiral" fields $ (as before) and $ = &e v obeying 
Va$ = and = (this last expression to be interpreted as a right 

action.). Evaluating the anticommutator yields: 

{V a , Va} = + Aite-^A^)) = iV a «. (42) 

The last equality defines the new object V a( j, whose lowest component is 
the usual covariant derivative. Taking one more commutator we obtain: 

W a = ~{V & ,Vaa\- (43) 

5 Not to be confused with the superpotential, also conventionally denoted with a W 
but without any spinor indices. 
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2.3 The complete Lagrangian 

Let us summarize all the results thus far by writing the most general J\f = 
1 gauge theory for a simple gauge group G. It is convenient to combine 
the usual Yang-Mills coupling g and the vacuum angle into a complex 
parameter 

t = — -i—. 44 
9 2vr 

Let be the matter chiral superfield transforming in a (generically re- 
ducible) representation of G. The complete Lagrangian takes the form 

S = — f d A xd 2 6 tr W a W a + — f d 4 xd 2 9 tr W & W & + 
16tt J 16vr J 

J d 4 xd 2 9d 2 $e v $ + J d 4 xd 2 8 W{$) + J d 4 xd 2 6 W($). (45) 

The superpotential W(3>) is invariant under G and all the coupling constants 
inside it are complex. In fact, in string theory it is quite natural to think of 
these couplings as the vacuum expectation values of some background chiral 
superfields as we will discuss later. 

It is instructive to write (|45[) in components in the WZ gauge. The gauge 
kinetic Lagrangian (first line in eq. ([45))) becomes 

J_ t r fA F ^ F ^ _ 2i\a tJ ~D fJi \ + D 2 ^ + JL e ^A trF/jiyjFpA . ( 46 ) 

The matter "kinetic" Lagrangian (first term in the second line of IJ45|) ) reads 

D^D^ - i^D^ + ff + V2i (0A^ - ^X4>) + 0D(f>, (47) 

where the indices 6 are contracted in the only possible way, i.e. (pXip = 
^) % \ aa T a \^) a j. Lastly, the contribution from the superpotential reads 

1 d 2 W , , dW . 

2 o4>id4)j o4>i 

Let us collect all terms contributing to the bosonic potential 

1 Q - dW dW - 

V = ~^D a2 - </>T a <i>D a - rfi ~ -Trrfi ~ j^f- (49) 
2g z a cpi d(p l 



6 The index a = 1 . . .dim(G) runs over the adjoint representation of G. The indices 
i, j = 1 . . . dim(R) label the representation R of $ and its conjugate representation R. 
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The auxiliary fields D and / appear only in these terms and we can solve 
for them as 



D a = -g 2 4>T a ^ fi 



dW 



f 



dW 



(50) 



dcpi 



Substituting the solution ijST))) in the potential gives 



dW dW 




(51) 



the last line to be interpreted as evaluated at the solution. 

The potential V is obviously non-negative, and the classical SUSY vacua 
are described by V = 0, i.e. by the simultaneous vanishing of the D and / 
terms, referred to as D-flatness and F-flatness conditions. 

The space of solutions to the D-flatness conditions cj}T a (j) = up to gauge 
transformations is an important object known as the classical moduli space 
of the theory. An extremely useful mathematical fact 7 is that such space 
can always be parameterized in terms of a set of independent holomorphic 
gauge invariants X r (rfi). One of the most common examples is a field in 
the adjoint representation of SU(N). In this case the D-flatness condition 
is cj)T a (/) = if abc (j) h (f) c = 0, that is, [<f> <fr] = as a matrix equation. The 
classical moduli space is parameterized by = tr (fi k for k = 2, . . . N. 

As for the F-flatness condition, it implies that when a superpotential 
is present classical SUSY vacua are found at its extrema. Thus, in the 
presence of a superpotential, some of the "moduli" X r defined above may 
be subjected to further restrictions. 

So far our discussion has been in terms of the classical theory and one 
must ask how this picture is modified in the quantum theory. We will see 
that, in general, the superpotential will be modified by (non-perturbative) 
quantum corrections but the property of the superpotential of identifying the 
SUSY vacua carries over to the low energy effective theory as long as one does 
not encounter any singularity, typically due to the presence of extra massless 
fields. In this review we will discuss at length how the superpotential is 
modified by quantum corrections but unfortunately we will not have time 
to discuss the many applications to the search for SUSY vacua. 

7 For a careful explanation accessible to physicists and a list of earlier references see |15|. 
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3 Basic Dynamical Facts 

Having seen the kinematic (classical) structure of SUSY gauge theories we 
now review some basic facts about their (quantum) dynamics. The aim of 
this section is to present enough information to be able to put the latest 
developments into context, without any ambition of being exhaustive. 

It is a familiar fact that in a quantum theory the couplings appear- 
ing in the physical quantities (observables) are actually running parameters 
depending on the renormalization scale [i and that they appear in the ob- 
servables in such a way that the total fi dependence cancels. All of this is 
still true in SUSY gauge theories but SUSY implies some simplifications in 
the running of the coupling constants. 

The simplest (and very important) example of this fact is the perturba- 
tive non-renormalization of the superpotential. This can already be seen in 
the WZ model, defined by the bare action 8 : 



Loop computations and renormalization can be done in a variety of equiv- 
alent ways - we use the so called "renormalized perturbation theory" and 
write action (|52[) in terms of the renormalized quantities: 



S = dVWZ$»+ dV 2 ^ Z ffl -$ 2 + Z A -$ 3 + c.c.= 



It is possible to develop Feynman rules directly in superspace without 
having to expand into component fields. Without presenting the details 
that are explained very clearly in the literature, the main point ^S] is that 
the most general term that can be generated in loop diagrams has only one 
Grassmann integral over all four 9: 



d 4 X! . . . d"x n d 2 9d 2 9G{ Xl . . . x^F^x^e, 9)... F n (x n , 9, 9), (54) 



where the Fi stand for products of superfields and their derivatives. We 
notice in passing that the expression would be valid if we also had the gauge 

8 We consider a cubic superpotential to ensure that the full theory is renormalizable in 
the ordinary sense. 




(52) 




(53) 
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superfield V with the possibility of the F{ being equal to V, D a V etc. in 

We know that the theory is renormalizable because it is so in components. 
Thus the only possible primitively ultraviolet (UV) divergent graphs are 
those with an external field structure already present in the bare Lagrangian. 
Suppose now that there were a UV divergent graph with only two (or three) 
external chiral fields <£>. The local counterterm needed to cancel it would be 
of the form (say in dimensional regularization, where e = 4 — d and I is the 
degree of the pole) 

j t J d A xd 2 6d 2 6 $ 2 (x,6,6) (55) 

but this diagram is identically zero since it is the integral of a chiral quantity 
over all of superspace. This means that diagrams with two or three chiral 
fields are always UV convergent and thus 5 m = 5\ = implying Z m = Z\ = 
1. At this stage, we do not need to worry about terms like 

r - D 2 

/ d 4 xd 2 9d 2 6 $, (56) 

which could be turned into a local F-term by replacing d 2 9 — ► D 2 and noting 
that D 2 D 2 & = — □<&, because such non-local behavior can never arise in the 
UV. 

The situation changes if we are interested in the effective action generat- 
ing the one particle irreducible diagrams (1PI action) ^7]. In this context, 
a UV convergent but IR divergent D-term of the type (|56j) can arise giving a 
term that looks like a finite contribution to the superpotential. One explicit 
example comes from the massless WZ model, where a diagram like Fig. Q 
gives a finite contribution to the 1PI effective superpotential: 

r \ \3\2 

Wxpi = J d 4 *d 2 fl + C(3)^^<I> 3 + . . . (57) 

and the dots indicate that there could be more contributions. 

Thus we see that the 1PI superpotential is (by definition) a holomorphic 
quantity in the fields but not necessarily in the coupling constants. This is 
an unpleasant state of affairs that suggests that we are not looking at the 
"right" quantity for our purposes. First of all, since we are interested in 
gauge theories, the 1PI action constructed via perturbation theory is not 
well defined since the fields appearing in the Lagrangian are not the correct 
degrees of freedom at low energy. Even more basically, if we think of the 
coupling constants as the vacuum expectation values (v.e.v.) of some very 
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Figure 1: The two-loop contribution to the one-particle effective superpo- 
tential. 



heavy chiral superfield ^H] ( as is natural in string theory), then a term like 
the last one in (|57|) is not supersymmetric and should be written instead as: 



where in the last line only the sum of all three terms is manifestly SUSY. In 
this framework, the whole expression should thus be thought of as a singular 
D-term. 

All this suggests that we forget about the 1PI action and look instead at 
the Wilsonian action defined as follows: If denotes the (bare) action at a 
scale no, then (for \i < /ip) denotes the action describing the same physics 
but where the degrees of freedom with momenta between [i and /xo have been 
integrated out. We can be more precise if we restrict our attention to the 
Wilsonian superpotential because any non-holomorphic dependence on the 
couplings must be regarded as a (now properly IR regularized) D-term |18j . 
The Wilsonian superpotential W e s is thus a holomorphic quantity in both 
fields and coupling constants and it is not perturbatively renormalized. Thus 
we write: 



where Wtree is the tree level superpotential and Wnon-pert a non-perturbative 
contribution that cannot be ruled out by the perturbative analysis above. 




(58) 



W eS = W tree + Wj 



non— pert > 



(59) 
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Such term is actually absent in the pure WZ model but will be present in 
many gauge theories and will indeed play a crucial role in all that follows. 

Let us make two more remarks. In the case where interacting mass- 
less fields are absent, the Wilsonian and the 1PI effective superpotentials 
coincide. Furthermore, the perturbative non-renormalization of the super- 
potential can also be proven using holomorphy and the symmetries of the 
theory as discussed in [T%] . 

Let us now turn to the renormalization of the gauge coupling constant 
r introduced in (|44|) . Since the gauge kinetic term is written as an inte- 
gral over chiral superspace one might be tempted to argue that the non- 
renormalization theorem would apply to this case too. However, this con- 
clusion is wrong because one can write 

J d 4 xd 2 6 tr W a W a = J d 4 xd 2 #d 2 tr (e~ v D a e v W a ). (60) 

The integrand is not gauge invariant, however the integral is (since it is just 
a rewriting of a manifestly gauge invariant one) and it is local. Thus, UV 
divergences of this type are allowed by the non-renormalization theorem and 
the coupling r is not protected from running. 

Having established that the coupling constant runs, we can use a piece 
of information well known from non-supersymmetric theories that states 
that the term, being the coefficient of a boundary term, does not get 
renormalized. Hence, the beta-function for r only involves Re(r) and thus 
cannot be holomorphic. The only exception would be if it were a constant, 
which would imply that g is renormalized only to one loop. Indeed, 

(61) 

where (3 is a real constant, implies 



d „ d 

M — 9 = and ^—g = —-^ g \ 62 
d\x djx 167r 



An explicit computation shows that this is not the case for the standard 
coupling constant that in general gets renormalized to all orders. Once 
again, this suggests that we should introduce a different coupling constant, 
to be referred to as the Wilsonian coupling, obeying (|61j). (|62|). The re- 
lation between the ordinary (1PI) and the Wilsonian coupling is of course 
non-holomorphic j 19 j and this explains why the two beta-functions disagree 
already at two loops. Even in this context, it can be shown that the lack of 
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holomorphy comes from an IR effect |19j and it is thus consistent to exclude 
it from the Wilsonian beta-function. From now on r will always denote the 
Wilsonian coupling. 

A standard computation then gives: 







^3r(Adjoint)-^T(p)^, (63) 



where T(p) denotes the index of the irreducible representations p in which 
$ transforms, normalized to T (Fundamental) = 1. 

Eq. (|61[) allows one to introduce a perturbatively RG invariant scale 
known as the holomorphic scale: 



,2 




A f3 = ^ e -2nr(p) = /^F^ 9 ( 64 ) 

Consider two theories related to each other by integrating out some field 
of mass m. Because the running is perturbatively only one-loop, the holo- 
morphic scales are easily matched as 



(65) 



where the tilde quantities refer to the IR theory where the field of mass m 
has been integrated out. 

The importance of the holomorphic scale A is that it controls the quan- 
tum corrections of the observables in the theory, with A — ► corresponding 
to the classical limit. 

Let us now return to the Wilsonian superpotential (|59|) . Notice that 
one can always write Wt Tee = ^ \ r X r ($i) for the gauge invariant quantities 
X r ($j) introduced in the discussion of the D-terms 9 . The question now is 
what is the form of VF non _ por t. One would expect that W n on-pert depends 
(holomorphically, of course) on the scale A, the couplings \ r and the in- 
variants X r ($i). But, in fact, W non _ per t is independent of the couplings A r . 
This fact is sometimes referred to as the "linearity principle" 20, because it 
implies that the full superpotential is linear in the couplings. 

The linearity principle is less surprising than it seems at first when one 
realizes that, having ruled out any perturbative dependence on A r , any non- 
perturbative dependence would have to have a singularity as A r — > in 



In some cases, some of the X r could also be dependent. 
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some direction in the complex plane. For instance, note that this rules out 
terms of the type e _1 / A that go to zero only if the origin is approached for 
i?e (A) > 0. The gauge coupling r is a different story because in this case the 
real part and the imaginary part enter differently (the real part is always 
positive) and thus the combination (|64|) appearing in A is allowed and it can 
be generated e.g. by instanton effects. 

Having established that W non -pcrt = ^non-pert (X r , A) it is sometimes 
possible to completely fix the functional dependence up to some numeri- 
cal constant by dimensional analysis and symmetry considerations. The 
most celebrated example of this is the Ameck-Dine-Seiberg superpoten- 
tial for SQCD, i.e. Supersymmetric SU(N C ) gauge theory with Nf quarks 
(1 < Nf < N c ) in the fundamental representation [21j . The new techniques 
which are the object of this review allow a more systematic derivation of 
these non-perturbative corrections to W e R applicable in more general situa- 
tions. 

What we have seen is that one can write 

W eS = AiXi + X 2 X 2 + ■■■ + W non ^ peit (Xx,X 2 , ... ,A U A 2 , ...), (66) 

where we have introduced more than one scale A s to allow for the case 
where the gauge group is not simple. Let us focus on, say, the dependence 
on Ai. Clearly, we can integrate out the field X\ at low energy (where the 
superpotential piece dominates) by solving 

wi w " = - (67) 

This is tantamount to doing a Legendre transform: 

d 

Al = — W n on-pert- (68) 

Solving (|68jl for X\ in terms of Ai and the other variables and substituting 
into (|66|) one obtains W e s where now the dependence on Ai will be com- 
plicated because we have integrated out its partner. Thus we can think 
of Ai and X\ as forming a canonical pair. One could integrate out all the 
composite fields X r to give W e g solely in terms of the couplings A r and the 
scales A s . 

The Legendre transform is clearly invertible and thus, contrary to our 
"Wilsonian" intuition, as long as we are only interested in F-terms we do 
not lose any information integrating out a field and, in fact, we can integrate 
it back in by reversing the procedure: 

{X r ) = ^-^eff- (69) 
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Note that this equation will be relevant later in a slightly different way: if 
we happen to know the v.e.v. of X r in some other way, eq. can be used 
to determine W e s as we will discuss in Section 7. 

So far the scales A s have played a passive role. However, the recent devel- 
opments depend crucially on the fact that we can also introduce a canonical 
conjugate for each A s . We can use the perturbatively exact definition of the 
holomorphic scale (|H4j) to rewrite (for each s) the gauge kinetic term as a 
(tree level plus one loop) contribution to the Wilsonian superpotential: 

^tr W a W a = (3 log(A//i)5, (70) 

l07T 

where S is the glueball superfield 

S = " 32^ tr ( W ° W ^ = 1^2 tr ( A ° A « + •••)• (71) 

The field S is a chiral superfield whose lowest component is proportional to 
the gluino bilinear tr X a X a . 

The linearity principle still applies in this case. While S does not appear 
in the original superpotential one can "integrate it in" by solving for A: 

<S> - fr±W«. (72) 

Notice that even without matter fields there is a superpotential for S: the 
so-called Veneziano-Yankielowicz (VY) superpotential [221 • Consider, as an 
example, pure SU (N) SYM theory for which (3 = 3N. By the assumption of 
confinement, we expect that all degrees of freedom are massive and thus the 
effective superpotential at low energy can only be a constant 10 . Dimensional 
analysis shows that W e g = aA 3 for some numerical constant a. Using l|72|) 
yields (S) = (a/N)A 3 . From an explicit instanton computation we know 
that 

(S) = A 3 , (73) 

that is, a = N. 11 

We can now express the non-perturbative superpotential as a function 
of S as 

^non-pert = ~NS log ^ + NS. (74) 

M 

10 A constant superpotential does not break rigid supersymmetry. 

11 The computation of the exact numerical coefficient is a subtle issue. To derive the 
correct result 17311 one needs to perform a computation at weak coupling. We refer to the 
extensive literature on the subject reviewed in \'23\. (See also I24|.') 
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Figure 2: The different energy scales 
Sometimes one also writes in a somewhat mixed notation 

W eff = ^ non _pe rt + SNS log ^=NS{l- log |jj = Wyy , (75) 

usually referred to as the VY superpotential. The advantage of this seem- 
ingly roundabout notation is that upon minimizing Wyy with respect to S 
one recovers ([73 [I. Thus in the following we will consider W e g as dependent 
on S although its "natural" variable should be A. 
The general situation is summarized in Tabled 
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Table 1: Fields and their corresponding couplings 



Given the F-terms as functions of one of the variables for each column, 
one can obtain the dependence on the other variables by performing a Leg- 
endre transform. 

Let us assume that the mass spectrum for these fields is as given in 
Fig. EJ Generically, all glueball fields S will be massive whereas we can 
split the remaining invariants into massive ones, schematically denoted by 
X' with couplings A', and massless ones X with couplings A. The "most 
Wilsonian" thing to do would be to express the F-terms as functions of the 
massless fields X, the couplings A' and the holomorphic scale A. 

What we will see in the following section however is that the new tech- 
niques allow us to determine the F-terms as functions of A, A' and S instead. 
Curiously, expressing the F-terms in this set of variables reveals simple de- 
pendences on other parameters of the theory such as the number of colors 
or the number of flavors that are obscured in the initial description. 
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4 Diagrammatic computation of the glueball su- 
perpotential 

In this section we will describe the first of the two main approaches to the 
calculation of the glueball superpotential, originally presented in [Jj- This 
approach is essentially perturbative in nature, and we will illustrate it by 
means of an example. We will follow the discussion of [Jj very closely and 
make some additional simplifications along the way for the sake of clarity. 

The computation of is based on the ideas put forward by Dijkgraaf 
and Vafa in pQ, where it was conjectured that the effective superpotential 
of a gauge theory with matter in the adjoint of U(N) and an arbitrary tree 
level superpotential Wt ree could be fully captured by a bosonic matrix model 
with a potential coinciding with Wtree- More specifically, one is instructed 
to compute the planar contribution 12 T to the free energy of the matrix 
model: 

e -^ T = [dMe-^ Wt ™( M \ (76) 



where M is a N x N hermitian matrix. Taking the large N limit while 
keeping S = eN constant singles out the planar diagrams and yields J- as 
a function of S and the couplings in Wt re e- The effective superpotential for 
the gauge theory 13 expressed in terms of these variables is then 

W eS = W VY + W D v, (77) 
where Wyy is the same as (|75|) and 

W DV = N-^F. (78) 

While the original conjecture was based on a chain of string dualities, the 
computation of is purely field theoretic, the reduction to a matrix model 
computation being a direct consequence of J\f = 1 SUSY in four dimensions, 
as we will show below. 

Consider the gauge theory described above and for simplicity take Wt r0 e 
to be a cubic superpotential: 

W tre e = ytr$ 2 + |tr$ 3 . (79) 

The reader should not confuse g with the gauge coupling that will always 
be "hidden" in A in the following. Note that renormalizability is not an 

12 We will see in detail how this is done. 

13 For simplicity we limit ourselves to discuss the unbroken phase at this stage. 
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issue here, since on the one hand the computation we will review is not 
sensitive to the degree of the superpotential, and on the other hand the 
theory we consider could be thought of as an effective field theory of an 
underlying renormalizable theory. Thus, the generalization to an arbitrary 
superpotential is straightforward. 

Thought as a polynomial in a complex variable z, the superpotential 
Wt V ee(z) = ^z 2 + |z 3 has two extrema 14 : z = and z = —m/g. The 
classical vacua are constructed by setting each (complex) eigenvalue of the 
N x N matrix (<3?) equal to one of the two extrema. This generically breaks 
the gauge symmetry from U(N) to U(p) x U(N — p). In this section we will 
study only the unbroken U(N) phase, i.e. we set the v.e.v. of $ identically 
to zero. 

The goal is to compute the effective superpotential W e s(S,g,m), as dis- 
cussed at the end of Section 3. Apart from the non-perturbative VY piece 
(|75j) already present in the pure gauge theory, the only contributions to W e g 
come from the perturbative evaluation of the path integral over the chiral 
superfield in a constant background for the vector superfield. 

Note that this is not in contradiction to the non-renormalization theo- 
rem, since when one integrates out S, all terms in W e g appear as powers of 
the holomorphic scale A. Thus it is true that W e R is a perturbative series 
in the coupling g, but each term is accompanied by a positive power of the 
essentially non-perturbative object A coming from the gauge sector. Were 
we to decouple the gauge sector from the chiral superfield, we would have 
to take A to zero, which would make W e a vanish. 

We choose the simplest possible background for our purposes, i.e. a 
constant gaugino X a and a vanishing gauge field, = 0. In order for the 
above background to satisfy trivially the equations of motion, we can set to 
zero the complex conjugate gaugino, X a = 0, employing the usual trick of 
thinking of X a and X a as independent objects. This is justified as long as 
we are interested in X a independent quantities such as the F-terms. 15 

In superfield terms, the background simply reads: 

V = -V2i9 2 eX, W a = -V2iX a . (80) 

14 Here and in the following we make an abuse of notation replacing $ by z in Wtree but 
refraining from taking the trace. 

It is possible to take a different approach and keep all generality in the background 
and in the superpotential, while using the full power of supergraph calculus (R.A. thanks 
D. Zanon for an interesting discussion on this issue). This would certainly be appropriate 
in an exhaustive review, however here we prefer not to introduce the more advanced 
material required for this treatment and instead we stick to a more intuitive approach as 
m El 
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There is a further restriction to be obeyed by the constant gluino background 
- a SUSY variation of the gauge field, SA^ = i^a^X, is required to be pure 
gauge, thus implying that the matrix valued \ a satisfy: 

{A a ,A /3 } = 0. (81) 

In this background, covariant spacetime derivatives reduce to ordinary 
partial derivatives, the gluino superfield W a is constant (and thus covariantly 
constant) and, most importantly, traces of more than two W a vanish. This 
comes about because eq. (|HTj) implies that, for instance, tvW a WpW^ must 
be totally antisymmetric in the spinorial indices, which take only two values. 
A further and less obvious consequence of (|8Tj) is that the glueball superfield 
S, as given in (|71|) . satisfies S h = 0, where h is the dual Coxeter number 
of the gauge group. This issue will be discussed in more detail in the next 
section, but it should be clear already at this stage that S is nilpotent 
because there are only a finite number of Grassmann variables in W a . 

We want to compute the holomorphic contribution to the partition func- 
tion: 

Z = j V<S>V$> e iS , (82) 

where the action is: 

S = j d 4 xd 2 0d 2 £tr (e- v $>e v $)+ J d 4 xd 2 9 WW$)+ J d 4 xd 2 9 W tTee (^)- 

(83) 

In the kinetic term we have written the adjoint action on $ explicitly, i.e. 
in terms of N x N matrices. Equivalently, one can define the adjoint action 
of V on <3? as Ad V<& = [V, <3?] . One can then exponentiate it by formal power 
series and check that: 

e AdV $ = l + Ady$ + ^(Ady) 2 $ + ... 

= l + [V,$] + ^[V,[V,$]] + --- = e v $e- v . (84) 

Both ways of writing this expression will be used. 

As far as the computation of the F-terms is concerned, we can set g = 
and m = 1 independently of g and m, since they are not going to enter 
any holomorphic quantity 16 . In this way the overall $ dependence will be 



16 Notice that setting fh = 1 will imply that some of the intermediate formulas that 
follow will not be manifestly dimensionally correct. The final (holomorphic) result will of 
course be. 
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quadratic. We can thus write: 

Z = j V<S> e' /d 4 ^ Wtree(#) [ £>$ e iS ° , (85) 

where 

S = J d 4 xd 2 0d 2 #tr (e- y $e v <S>) + J d A xd 2 6 ^tr$ 2 . (86) 

We start by evaluating the (Gaussian) path integral over 
The action So can be reformulated by defining: 

$ = e -v$ e v = e -Adv^ = ^Adv (8?) 

Using the covariant fermionic derivatives defined in (|41jl , where V should be 
intended in the adjoint representation (i.e. V a = e~ AdV D a e AdV ) the fields 
now obey Va<& = and = 0. More importantly, the path integral over 
<& can be safely replaced by a path integral over <5, since no anomalies can 
arise, <3? being in the adjoint representation. 

Finally, to complete the square, we need to turn the integral over an- 
tichiral superspace into an integral over full superspace: 

So = J d 4 xd 2 9d 2 6 tr ($<5>) + J d A xd 2 6 ^tr$ 2 

d A xd 2 9d 2 6 tr ^<M» + ^$^2$ 

where we have used the fact that J d 2 can be replaced by V 2 when inside 
the spacetime integral and the trace. 

Note that V 2 is not invertible in all of superspace (for instance it vanishes 
on covariantly antichiral superfields like $), but we can define its inverse 
up to an element of its kernel. In the literature, it is also often written 
equivalently either as (\7 2 \7 2 )^ 1 \7 2 or as V 2 (V 2 V 2 ) -1 . The first expression 
is obviously the left-inverse of V 2 , while the second is the right-inverse, and 
again they coincide up to an element of the kernel as can be checked by 
conjugation. 

Completing the square, the action can thus be rewritten as: 



So = I d A xd 2 6d 2 9 x - 



tr ($ + V 2 $)^r($ + V 2 $) - tr (3>V 2 $) 



(89) 
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Now the first piece is trivially a quadratic integral on antichiral superspace 
for the shifted antichiral superfield = $ + V 2< £, while the second piece 
can be turned into an integral over chiral superspace: 

S = j d 4 xd 2 itr ($') 2 " J itr ($V 2 V 2 $). (90) 

The path integral over <3? can be harmlessly turned into a path integral over 
which is Gaussian and independent of the background, thus contributing 
just a constant to the partition function Z. 

Returning to Q85|). we are left with a chiral path integral over <£. Con- 
trary to the integral over the action here contains interactions, forcing 
us to a perturbative evaluation. We will however encounter a remarkable 
simplification due to supersymmetry which will lead us to consider a matrix 
model. 

To obtain the explicit form of the propagator for <3? we need to expand the 
expression V 2 V 2 <1?. Recalling the definitions, we have that (in this simple 
background): 

V 2 V 2 $ = D 2 e- AdV D 2 e AdV <S> = -□<£> - AdW a D a $ = -□<£> - {W a ,D a $}, 

(91) 

where in this case the adjoint action gives rise to an anticommutator. 

The path integral (|85|) is thus reduced to the following chiral path inte- 
gral: 

Z' = j V$ e* f d4xd2e |tr$(n+AdH/ Q D Q +m)$+ftr<I> 3 ^ g2 ) 

The quantity Z' defined above is the purely holomorphic part of the full 
partition function (|H2|). i.e. the piece that depends only on the holomorphic 
couplings and on the holomorphic background X a . We can thus associate it 
with the coupling dependent part of the effective superpotential: 

Z> = e iJ d4 * d2e W DV (S,g,m) _ (-93) 

The procedure outlined above is quite general, and it applies to any 
group, matter content, and tree level superpotential. 

There are however two important points to make. The first one is that, 
since Z' is the partition function obtained by integrating over the matter 
fields, there is no hope of getting the piece of the effective superpoten- 
tial pertaining to the pure gauge low-energy dynamics, i.e. the Veneziano- 
Yankielowicz piece (|75|l. in this way. This piece is independent of the matter 
couplings g and m and cannot arise from a perturbative evaluation of (|92j). 
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The second point is that the perturbative evaluation of ()92|) has to take 
into account the fact that the classical gluino background satisfies (|81j) . 
As said before, this implies that the gauge invariant gluino bilinear S is 
nilpotent, i.e. S h = 0. For instance, in the U(N) case considered above, 
Wpy (S, g, m) would have a truncated Taylor expansion which stops at S N . 
Of course, this makes no sense when S is considered as an effective chiral 
superfield, and we will indeed find superpotentials Wdv(S, g,rn) with an 
infinite Taylor expansion, of which the first h — 1 terms coincide with the 
ones computed by the method above. 

In some cases however, one can compute Wdv{S, g, m) to arbitrary pow- 
ers of S without any obstacle. This is possible for instance in the case con- 
sidered here because the dependence on N factorizes in front of Wdv, as 
stated in (|78|). Since T is independent of N we can compute its power se- 
ries expansion to arbitrary order by formally considering N large and then 
substitute into (|78|) with the actual value of N. In other cases discrepancies 
with exact field theory results can show up at order h, as discussed in [2*5] . 

We now proceed to evaluate perturbatively (|§2"|) and to show that in the 
case of U(N) with adjoint matter, it reduces to the evaluation of a matrix 
model integral. 

First of all, it will be useful to go to momentum space. Following [Jj, we 
introduce not only the usual (bosonic) 4- momentum p^, but also a fermionic 
momentum ir a conjugate to 6 a . In the quadratic part of the action in (|92[). 
we Wick rotate to Euclidean space and Fourier transform in both the bosonic 
and fermionic coordinates by substituting: 

□ -► -a Eud -> p 2 , D a -» TT a . (94) 

It may be surprising that we do not include a factor of i in the last definition 
but this is consistent with our choice of hermiticity for D a . 

We now turn to the propagator (<I>"(p, 7r)<3?^(— p, — 7r)) for <F In the 
following $ is thought of as a N x N matrix and the roman indices a,b, . . . 
run from 1 to N. Looking at the kinetic term in (j92[) we read off the inverse 
propagator: 

F? d (p, vr) = ( P 2 + m)5 a d 5 b c + (W aa d 5 b c - W ab c 5^ a . (95) 
We can write the propagator itself using a Schwinger parameter s: 

A(p,7T)= / ds e s(p 2 +m+AdW^ a ) ^ ( 96) 

Jo 
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Figure 3: The parts Yf d and G^ d of the propagator and the cubic vertex. 



where, as usual, the exponential is defined by power series expansion 17 . 
Actually, since n a is a Grassmann two-component spinor, the expansion 
stops at the second order. Thus we have: 

e -«G = 1 _ fl G+-* 2 G 2 , ( 9? ) 

with 

lf d ee 6 a d 5 b c , Gf d = (W aa d 5 b - W ab c 5 a d )* a . (98) 

In order to compute Z' we have to evaluate vacuum diagrams of the 
theory defined by the action in (|92j) . Because the matter field <I> is in the 
adjoint, the diagrams can be written using the double line notation |26| 
(see Fig. OJ). The identity in (|97j) is represented by two lines with opposite 
orientations, while the G and G 2 terms are represented by the same lines 
with one or two W a insertions respectively. The cubic vertex, being a trace, 
is represented by three lines joining three different two-line propagators. 

Having established the above, one can assign a topology to any given 
graph in the usual way. The graph represents a discretization of a two- 
dimensional closed surface over which it can be drawn without any twisting 
or intersection. Call V the number of vertices, P the number of propagators 
and F the number of faces i.e. regions encircled by a line, including the 
external one. Then the Euler characteristic of the surface is given by: 

x = V - P + F, (99) 

and it is a topological invariant. For instance, a sphere has Euler charac- 
teristic x = 2 and a torus has x = 0- Actually, every handle one "adds" to 
the surface reduces x by 2, so that the following formula is valid for closed 
orientable surfaces: 

X = 2 - 2h, (100) 
"Notice that {W a ,n a ®} = \W a ir a ,$\. 
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where h here is the number of handles, also referred to as the genus of the 
surface. 

Allowing for open and/or unoriented surfaces, the Euler characteristic 
can take odd values (but always keeping \ < 2 for a single surface). For 
instance, \ = 1 for a disk (open, oriented surface) and for the projective 
plane (closed, unoriented surface). 

An important result that is obtained even before starting any compu- 
tation is that only planar graphs contribute to the superpotential. This can 
be understood as follows. A generic (connected) vacuum diagram has V 
vertices and P propagators. Every vertex carries with it an integral over 
chiral (coordinate) superspace, while every propagator, after Fourier trans- 
forming it, carries an integral over chiral momentum superspace. All but 
one of the coordinate superspace integrals give delta functions enforcing 
momentum conservation (both bosonic and fermionic) at every vertex, the 
remaining overall chiral superspace integral being the one which appears 
in the effective action in ()93j) . Solving for the delta functions, one is left 
with P — V + 1 = L momentum superspace (loop) integrals. So far this is 
just ordinary loop counting applied also to the Grassmann momenta. Thus 
one has to integrate over 2L Grassmannian momenta ir a . In the integrand, 
these momenta necessarily appear in bilinears together with the gluino back- 
ground W a because of the specific form of the propagator $Fg§ . Since traces 
of more than two W a vanish due to (|81[) , the result of the integration is non- 
zero only if we are able to write at least L traces. Remembering now that 
the number of traces coincides with the number of closed lines, or "faces" 
of the graph, we get the following constraint on the graph's topology: 

F>L = P-V + 1 x>l- ( 101 ) 

This means that in the case of the matter field in the adjoint of U(N), 
we are concerned only with graphs with the topology of a sphere, that is 
planar graphs. All the other ones are exactly zero already at finite N. In 
other words, planarity is an exact consequence of supersymmetry, and not 
a leading order approximation due to a large N limit. 

When open or unoriented graphs have to be taken into account, also 
graphs with one boundary or one cross cap give a non- vanishing result. 
Graphs with boundaries arise in theories with matter fields in the funda- 
mental representation, whose propagator has only one color line (possibly 
supplemented by a flavor line), while unoriented graphs arise in theories 
where the gauge group is SO(N) or Sp(N). Note also that graphs of higher 
genera become important when the theory is coupled to gravity, in the back- 
ground of which they no longer vanish. 
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Let us now compute in detail one diagram in order to give a feeling of 
how a field theory computation can be reduced to a matrix model one. 
Let us compute the "stop sign" diagram of Fig. |IJ 



Figure 4: The stop sign diagram. 



This diagram yields: 

y y"^^d 2 7r 1 d 2 7r 2 A^(p 1 , 7 r 1 )A^(p 2 , 7r 2 )A^(- Pl - p 2 , -tq - vr 2 ) 

(102) 

2 , 

The overall factor of is understood as follows. We have a from the Taylor 
expansion to second order of the exponential of the interaction, a {g/3) 2 from 
the two interaction vertices, and a 3 from the number of independent ways 
one can contract the <3? to give a planar graph 18 . 

In terms of the three Schwinger parameters, we can rewrite the diagram 

as: 

^ J d Sl ds 2 ds 3 e- m ( Sl+S2+ ^/ bose / fermi . (103) 

The bosonic integral is straightforward, since it does not imply any group 
theory factors. It is given by: 

_ /d^d^ 2_ S2p 2_ S3(pi+p2)2 

Vbosc - J (2?r)4 (27r)4 

' ' (104) 



(47T) 4 (sis 2 + S1S3 + s 2 s 3 ) 2 ' 
The computation of the fermionic integral is more involved, since the inte- 



18 There are also three more contractions that do not yield a planar graph and thus 
vanish according to the above argument. This explains the discrepancy between 11021 
and the factor that would arise in a single component 3 theory. 
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grand involves traces over the gauge group: 

Jfenni = [d\id\ 2 (e-'^^y(e-^ w ^)^(e^ wa ^^Y f 

J V J cd\ ) af\ J eb 

(105) 

The integral over the fermionic momenta is saturated only if we bring down 
two 7Ti and two 7^: 

J d 2 Kid 2 -n 2 A.A = 1. (106) 

The result thus gives an overall fourth power of the gluino W a . These 
four gluinos have to be inserted at some point of the three color loops of 
the diagram. Since we cannot put more than two gluinos on the same line 
(otherwise the trace would vanish), the two choices we are left with are 
either to put two W a on two color loops and leave the third loop without 
insertion, or to put two W a on one loop and one each in the two remaining 
color loops. 

The latter combination leads to a term including tr VF a tr W a , which is 
the effective kinetic term for the decoupled U(l). Terms of this type are 
interesting in the case of broken gauge group but here we focus only on the 
effective superpotential for the gluino bilinear S, i.e. we are concerned only 
with those combinations where either two W a or none are inserted in each 
color loop. The U(l) couplings will be briefly discussed in Section 7. 

It is easy to work out all the insertions pictorially, instead of developing 
the propagators by brute force. One can then realize that two different kinds 
of diagrams contribute. The first, depicted on the left of Fig. [SJ arises when 




Figure 5: The stop sign with two W a insertions at two of the propagators 
(left) and the stop sign with one W a insertion at two of the propagators and 
one insertion of two W a at the third propagator (right). 

two propagators have two insertions and the third has none 19 . Diagrams of 
19 At this point it is important not to confuse propagators and color loops - propagators 
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these type contribute to Ifermi a quantity: 

3iV(tr W«W a ) 2 (±s 2 \s 2 + \s\\sl + \sl\s^ , (107) 

where the factor of N is given by the color loop without insertions, tr 1 = N, 
while the factor of 3 comes from the choice of this color loop. 

The second kind of diagram, depicted on the right of Fig. arises 
when one propagator carries two insertions and the other two one each. The 
contribution from these diagrams is: 

3JV(tr W a W a f Qs 2 s 2 s 3 + si~s 2 2 s 3 + s^s 2 ^ , (108) 

the numerical prefactor having the same origin as above. 
Summing these contributions, we get, using (|71|): 

/fcrmi = 3iV(47r) 4 5 2 ( S is 2 + Sl s 3 + s 2 s 3 ) 2 . (109) 

Comparing now with the result of the bosonic integral (|1U4|) , we observe the 
striking fact that the Si dependent part of the numerator and the denomi- 
nator exactly cancel 20 , leaving: 

-fbose-ffermi = 3N S 2 (HO) 

that is, a result independent of the Schwinger parameters Si ! 

Plugging this back into ()103|) . the integral over s\, s 2 and S3 becomes 
trivial: 

[ d Sl ds 2 ds 3 e - m (si+S2+s 3 ) = _L (in) 

J m A 

In other words, all that is left from the propagators is a contribution of -jr 
each, exactly what one would expect from a zero-dimensional field theory, 
i.e. from a matrix model. In particular, note that the result is finite. This 
is non-trivial since after all we are dealing with a two-loop computation. 
The final result for the stop sign diagram is thus: 

1 a 2 

-■^3NS 2 . (112) 

6 m 6 

To order g 2 there is only another diagram, shown in Fig. El contributing 
to the effective superpotential. We leave it as an exercise to check that, after 

can carry at most two insertions because their expansion stops at second order. 

20 Curiously, this works only in four spacetime dimensions, as the overall powers of Si in 
the numerator are determined by the dimension of the Weyl spinor, whereas the ones of 
the denominator are determined by the dimension of spacetime. 
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Figure 6: The drum-bell diagram. 



all the W a insertions are made, it contributes a quantity 

1 a 2 

-^3NS 2 . (113) 

Summing (|112j) and (j!13[) and recalling the overall minus sign coming from 
the Wick rotation we get the full contribution to the effective superpotential 
to this order: 

a 2 

W DV = -2N^S 2 + ... (114) 
mr 

In all generality, we can argue as above that for a planar diagram with 
F color loops, the S = S F ~ 1 term will have a prefactor of N from the color 
loop without insertions and a factor of F from the choice of such a loop: 

4osc/fermi = FNS^ 1 = N^S F . (115) 

The numerical and coupling dependent prefactor will depend on the inter- 
action vertex and on the combinatorics. That the Si dependence completely 
cancels in the above expression is less trivial to see, but can actually be 
proven in all generality introducing a new set of 2L auxiliary Grassmann 
coordinates to enforce the fact that one allows for at most two insertions of 
W a in any color loop. We preferred here to present a sample computation, 
and we refer the interested reader to the original derivation [7]. 

The relevant question now is whether we can define a matrix model 
from which to extract all the information we need to compute the effec- 
tive superpotential of the gauge theory. The answer turns out to be quite 
straightforward. 

Let us look at the following matrix model: 

Z = jdMe-~^^ M2+ i Mi ). (116) 
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Here M is an N x N hermitian matrix, where N is not related to the N of 
the gauge theory. Actually, for the time being the only explicit link between 
the matrix model and the gauge theory is that the potential of the matrix 
model is given by the tree level superpotential Wt ree of the gauge theory, 
after substituting <I> for M. Otherwise, note for instance that the matrix 
model is purely bosonic. The parameter e is introduced for later convenience, 
and plays the role of an overall coupling constant. 

Strictly speaking the potential in (|116|) is not well defined because it 
does not have a minimum. This should not be however a matter of worry, 
since we will use (j!16j) to define a (formal) perturbative expansion around 
its local minimum (recall that we are interested in the case where, on the 
gauge theory side, the v.e.v. of $ does not break the gauge group). 

Indeed, we write: 




(117) 



In the second line above, we have used the correlation functions of the "free" 
matrix model, that is the one given by a purely quadratic potential. The 
free propagator is given by: 

{M£M c d ) = U a d 5l (118) 

and can also be represented by a double line as in the previous set up. Note 
however that here the roman indices run between 1 and N. 

The vacuum diagrams of the interacting theory are represented by the 
correlation functions (|117|) of the free theory. Every vertex is multiplied by 
a factor of g/e. Hence a generic diagram will give a contribution of: 




(119) 



where c is a numerical coefficient including the powers of | coming from the 
vertices, the coefficients of the Taylor expansion of the exponential and the 
various combinatorics. Again, F is the number of index loops, i.e. traces. 

The matrix model simplifies greatly when one takes the large N limit. In 
order to do this properly, let us take simultaneously the limit e — > 0, while 
keeping the combination: 

eN = S (120) 
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constant (it is of course not an accident that the constant is called S, we 
will shortly see below why). 

We now solve for N in (|119|) . and recalling the expression for the Euler 
characteristic (j99|) . we obtain for the graph: 



c e 



-* (m~ p g v S F ) . (121) 



The expression between parenthesis is finite in the limit, and when e — > 
the planar diagrams (x = 2) will clearly dominate, the other ones being 
suppressed by powers of 1/N 2 . 

Reexpressing the partition function for the matrix model as: 

Z = e -^X=2-^"x=0-^x=-2 + - ) ( 122 ) 

we observe that the result of computing all planar diagrams is contained in 
the planar free energy T = J- X= 2(S, g, m). For instance, the contribution to 
T from the "stop sign" diagram considered before, is: 

-l^s9 2 S 3 , (123) 

where the power of S is now simply the number of color index loops. The 
cubic matrix model can be solved exactly in the planar limit, as first shown 
in jlUj and in Appendix B we indeed show how to obtain the complete T . 

Recalling the expression (|115|) derived on the gauge theory side, the only 
non-trivial pieces of informations to compute are the numerical factors in 
front of every term in the perturbative expansion in S. These are precisely 
the numbers contained in the function T '■ Hence the effective superpotential 
of the gauge theory is given in terms of the matrix model planar free energy 
by: 

W DV (S, g, m) = N-^HS, 9, m). (124) 

Let us underscore again the important point - planarity is achieved through 
the large N limit on the matrix model side, but is an exact feature at 
finite N on the field theory side. The reason is that on the field theory 
side non-planar diagrams necessarily include traces of more than two gluino 
super fields, and these vanish exactly. 

Furthermore, the only way N enters in the effective superpotential is 
through a multiplicative factor in (|124j) . Thus, we see that, even if on the 
gauge theory side we should trust the effective superpotential only up to the 
term S N , the matrix model instructs us to go on and compute any diagram 
in the same way. This works in this case because the functional form of 



36 



the superpotential does not change when N is varied, and thus one could 
also take a large value of N on the gauge theory side. This generalization 
must of course rely on considering S a dynamical effective field (the glueball 
superfield). 

We have already stressed that the perturbative approach does not give 
us a systematic way to compute the pure gauge part of the effective su- 
perpotential, namely the Veneziano-Yankielowicz superpotential for a pure 
SU (N) SUSY gauge theory. This must be added to the DV term in order 
to obtain the full effective superpotential. The total effective superpotential 
thus reads: 

W eS = W VY + W DV = Ns(l- log j^j + N^(S, g, m). (125) 

The holomorphic scale A appearing in (|125|) is the one for the low-energy 
pure SU(N) theory, when the massive matter field $ has already been in- 
tegrated out. Using the matching of scales (|65[). we see that it is related to 
the high energy scale A by A 3 = mA 2 . Were we to be strict and use only the 
high energy scale in the VY superpotential, the matching of scales would 
come about from the vacuum diagram of the free matrix model, which is 
given by S 2 logm. 

It may sound a bit unpleasant that one must add the VY piece "by 
hand". However, a very interesting fact well emphasized already in the 
earlier literature EZl is that the VY superpotential can be fully recovered 
in the matrix model if one takes into account the volume of the U(N) group 
rotating the hermitian matrix M. Essentially, one can show that: 

(1) oc (Vol(U(N))) ~* ~ log^+-, (126) 

so that after substituting for S, the leading term looks like a contribution 
to the planar free energy, and then the typical S log S term also arises from 
(|124|) 21 . Although this observation can be justified by string duality it is 
difficult to find a field theoretical argument for it. In the framework we have 
been working, we have consistently taken (1) = 1 and added the VY piece 
by hand. 

21 After the first version of this review appeared, the matrix model prescription which 
yields exactly the VY superpotential was spelled out in |28|. 
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5 The chiral ring 



One of the uses of the effective superpotential is that of determining the 
vacuum structure of the theory, namely the expectation values of certain 
gauge invariant operators. It thus comes natural to define an equivalence 
class of such operators identifying those that have the same vacuum expec- 
tation value in a SUSY vacuum. This leads to the concept of the chiral ring 
that we now turn to discuss. 

What we have seen in Section 2 is that a chiral superfield is defined by 
D a <& = and its lowest component 4>{x) is annihilated by Q a : 

[Qa,<f>(x))=0. (127) 

The same facts are true for the gluino X a (x) which is the lowest com- 
ponent of the chiral superfield W a , only now we must switch commutators 
and anticommutators to account for the fact that X a (x) is a fermion: 

{Qa,K(x)} = 0. (128) 

We call chiral operators those gauge invariant operators 0{x) that are 
annihilated by the Q a . For, say, a bosonic operator: 

[Q a ,O(x)] = 0. (129) 

As we have seen before, given such an operator one can always generate 
its superpartners by acting on it with the Q a to fill up a chiral multiplet. 
Vice-versa, given a gauge invariant chiral superfield its lowest component 
will be a chiral operator. 

Before going on studying the properties of such operators, we should dis- 
cuss how to construct them in a gauge theory. Here we have at our disposal 
the gauge variant objects (p(x) and X a (x) that are annihilated by Q a . One 
could thus construct gauge invariant composite operators, such as tr X 2 (x), 
that are naively annihilated by Q a by using the Leibniz rule. But one must 
be careful when defining composite operators in quantum field theory. In 
particular, when using point-splitting regularization, one must introduce a 
Wilson line to preserve gauge invariance, as we will discuss shortly. One 
then has to check that such regularization does not spoil eq. (|129|) . In this 
case it turns out that the naive expectation is correct. 

Consider, for concreteness, two fields q and q, the lowest components 
of the quark chiral superfields in SQCD, transforming in the fundamental 
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and antifundamental representations of the gauge group. Define the gauge 
invariant composite operator 



qq{x 



) 



lim q(x + e/2)P exp -i / A^(x)dx^ q(x - e/2) 

V Jx-e/2 J 



lim g(x + e/2) (1 - ieM^x) +...)q(x- e/2), 



(130) 



where the explicit insertion of the Wilson line guarantees gauge invariance. 
The dots in the second expression represent the high order terms in the 
Wilson line needed to preserve gauge invariance to all orders but they will 
not be needed here. 

When acting with on the LHS of ()13U|) we get zero by definition when 
hitting q and q but the SUSY variation of gives rise to a gluino (cfr. eqs. 
(gj) and (O): 



[Qa,qq(x)] =—~ lim q(x + e/2)e^-X a (x)q(x - e/2). (131) 



The quantity in the limit is naively linear in e but we could get into trouble 
if the OPE between q and q had a singularity in e so that a non-zero result 
could arise. Fortunately, no such singularity is present in this case, as can 
be seen by noticing, for instance, that the propagator between q and q is 



It is a general fact that there are no singularities in the OPE of lowest 
components of chiral fields and thus the gauge invariant composite operators 
made out of them obey the naive equation [Qa,0(x)] = and are thus 
chiral operators. We will see later, when discussing the Konishi anomaly, 
that there are other cases (not just involving 4>{x) and \ a {x)) where such 
singularity arises and the naive SUSY transformations are modified. 

Now that we have constructed chiral operators, we can immediately de- 
rive the important consequence jSHj that the v.e.v. of an arbitrary (time 
ordered) product of such operators is totally independent of their spacetime 
position. Consider for instance the product of two bosonic chiral operators 
0%(x) and 02 (y) and take the derivative with respect to x^: 





(132) 




d 



O 1 (x)O 2 (y))|0)+^ 



(0\[O 1 (x),O 2 (y)]\0)5(x°-y ). 



dx^ 
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Both terms vanish separately. The first because (for, say, x° > y°): 

<O|^0i(aOCMi/)|O) = -i(0\[P l ,,O 1 (x)]O 2 (y)\0) (133) 

= l -^ a (0\[{Q Q ,Qa},O 1 (x)]O 2 (y)\0) 
= l -a^(0\{Q £e ,[Q a ,O 1 (x)]}O 2 (y)\0) 
= ^(0\{Qa,[Qa,Oi(x)]O 2 (y)}\0) 
= 0. 

The idea behind the manipulations above is to use the SUSY algebra, the 
Jacobi identity and the chirality of the operators to bring Q„ to act on the 
vacuum state, where it vanishes assuming that the vacuum is supersymmet- 
ric. 

The second term in (|132|) is also zero because the equal time commutator 
vanishes since the same arguments leading to (|133(l apply more generally to 
the OPE of the two operators. 

Given that the v.e.v. of chiral operators is independent of their position, 
one can go to the limit of large separation and apply cluster decomposition 
to obtain: 

<O|r(0i(a;i) . . . O n {x n )) |0> = (d) . . . (O n ), (134) 

where there is no longer any need to specify the positions. This is the 
important property of factorization of correlation functions involving chiral 
operators. 

Since objects of the type {Qa, ■■■} do not contribute to the expectation 
values in a SUSY vacuum, it is natural to define an equivalence relation be- 
tween chiral operators. Two chiral operators 0\{x) and 2 {x) are equivalent 
if there exist a gauge invariant operator Xa(x) such that 

01 (z) = 2 (x) + {Q ( *,X A (x)}, (135) 

where the generalization to operators carrying Lorentz indices should be 
obvious - the extra indices just appear in the same way in Oi(x), 2 {x) and 
Xa{x) and the anticommutator in (|135l) becomes a commutator if 0\{x) 
and 2 (x) are fermionic. 

Note that for (|135|) to be consistent, Xq(i) must satisfy: 

[Qg,{Q™iXa(x)}} =0 (136) 



40 



which is not automatic from the SUSY algebra. 

The set of equivalence classes of chiral operators under ()135j) forms a 
ring, known as the chiral ring. It is easy to check that the product of equiv- 
alence classes is well defined (i.e. it is independent of the representative) by 
applying the Leibniz rule a few times. 

The equivalence (j!35j) can also be formulated in superspace by saying 
that, if 0\{x) and 02 0*0 are the lowest components of a chiral superfield Si 
and £2 respectively, there exist a (gauge invariant) superfield Y& such that: 

Sx = S 2 + B*Y & . (137) 

We will also refer to two such superfields as equivalent in the chiral ring. 
Again, the superfield Y& must obey D^D a Ya = for consistency. Fur- 
thermore, the relation D&Yg + D^Y^ = can always be imposed since it 
disappears from ()137|) . These two conditions together allow one to "solve" 
the constraint as Y& = D&Z, for some superfield Z in terms of which (|137|) 
now reads 

Si = £ 2 + D & DaZ = £ 2 + 2D 2 Z. (138) 

Thus, two superfields are equivalent if they differ by a "chirally exact" term. 
We have also shown, in passing, that the lowest component of a chirally exact 
object is the SUSY variation of some other operator and, thus, its v.e.v. is 
identically zero in a SUSY vacuum. 

Let us make all of these statements concrete and apply the above for- 
malism to gauge theories. Let <3? be any chiral superfield transforming under 
the representation T a of the gauge group. We write W a = W^T a . As usual, 
we denote the lowest components of $ and W a (in the WZ gauge) by (f) and 
—\^2i\ a respectively. The basic relation that we will extensively use is: 

A^ = -^=[Q d ,V ad ^], (139) 

or, equivalently, in terms of superfields 

W a <5> = V 2 V Q $. (140) 

Since the fields appearing in the Lagrangian are gauge variant we must use 
covariant derivatives. Of course, when acting on gauge invariant combina- 
tions, V 2 becomes D 2 . 

The proof of these relations is simple: e.g. 

W a § = ^[V A , V ad ]$ = l -V«V a ^ = Jv d {V d , V Q }$ = V 2 V Q $. (141) 
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(In the case where $ is in the Adjoint representation, it is customary to 
think of it as a matrix and write [W a , $] instead.) 

We now use (|140j) to obtain relations in the chiral ring. We call these 
relations "classical" , and similarly for the ring generated by them, since they 
receive quantum corrections as we will discuss later on. 

Let us specialize, for concreteness, to pure gauge theories, where the only 
chiral superfield is W a itself, transforming in the adjoint representation. Eq. 
(fTlUl) reads: 

{W a ,Wp} = V 2 V a Wp. (142) 

Note that the RHS of Q142|) is symmetric in a and (3 because of the relation 
1 \J a W a = V^PF", which can be easily checked using the definitions of W a 
and W a . 

Consider now taking the trace of a string of W a . Because of (|142[) . 
anticommuting any two W a will contribute only a trivial element in the 
chiral ring, e.g.: 

tr (WaWpWj) = -tr (WaWyWp) + D 2 tr {W a VpW y ). (143) 

Using again l|142|) , we see that the only non-trivial piece of the trace is when 
all the indices are antisymmetrized, thus excluding all traces of three or 
more W a . 

Hence, of these type of invariants, the only one that is non-trivial in the 
chiral ring is the glueball superfield introduced in (|71|): 

S = 3_tr W a W a . (144) 

Notice that, in spite of the fact that S oc D 2 ti (W a e~ v (D a e v )), the glueball 
superfield is not chirally exact because D a tr (W a e~ v (D a e v )) is not gauge 
invariant. 

Let us now be more exhaustive and discuss the full chiral ring of all pure 
gauge theories based on simple Lie groups. 

For a pure SU (N) gauge theory, eq. IJ144J) is the end of the story because 
traces of W a exhaust all the invariants one can construct. In fact, the only 
primitive invariants at one's disposal to make a singlet are Sj, €i lt ...i N and 
e iIr " !jv , but €j!,...j N x e* 1 '— *jv [ s proportional to products of S) and one is only 
left with traces. 

For other gauge groups one has to be more careful. For instance, for 
the remaining classical groups one could use the extra primitive invariants 
((P- 7 = 5 3% for the orthogonal and r\ %:! = —rf 1 for symplectic groups) to raise 
one index of the W a and then try to use one €i ± i N to make a singlet. But 
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even this will not work: For SO{2k + 1) the e-tensor has an odd number 
of indices whereas one needs an even one. For Sp(k) Wa is symmetric 
and thus vanishes when contracted with e and, finally, for SO (2k) (k > 2) 
the contraction of two W\ or of two W2 with e vanishes because of their 
Grassmann nature. 

Thus, for pure gauge theories based on classical groups the chiral ring is 
generated by the glueball superfield [3111 • This is conjectured to be true even 
for exceptional groups although the presence of further primitive invariants 
(e.g. a l i k for G2) makes this statement non-trivial. Amongst the exceptional 
groups, so far the conjecture has been proven only for G2 [3*T] . 

Having established the generators of the chiral ring, we now proceed to 
see that it has a further structure. An important and non-trivial fact that 
one can show using group theory is that relation (|142|) also implies 

S h = in the classical chiral ring. (145) 

The integer h is the dual Coxeter number of the group, i.e. half of the 
index of the adjoint representation. Eq. (|145|) has been proven for the 
classical groups (HOI and for G2 |S] and it is believed to be true for the 
other exceptional groups as well. The list of dual Coxeter numbers is given 
by Table 03 We will not go through the proof of this fact in general but just 



SU{k) 


SO(2k) 


Sp(k) 


SO(2k + l) 


G 2 


Fa 


Ee 


E-; 


Es 


k 


2k -2 


k+1 


2k -1 


4 


9 


12 


18 


30 



Table 2: Different groups and their dual Coxeter numbers. 



notice that for SU(2) it can be easily obtained from the relation (valid for 
any four SU(2) generators) 

tr (ABCD) = - (tr [AB)tv (CD) + tr (AD)tr (BC) - tr (AC)tr (£>£)) . 

(146) 

Substituting A = B = W\ and C = D = W2 and using Fermi statistics one 
gets 

S 2 oc (tr(WiW 2 )) =tr(W 1 W 1 W 2 W2), (147) 

the last term being zero in the classical chiral ring. 

What we have obtained so far is very similar to what was obtained in 
the previous section where imposing the condition 1)81(1 on the background 
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has set to zero precisely all the operators that now are shown to be trivial 
in the classical chiral ring. 

However, if eq. Q145|) were correct at the full quantum level, it would 
imply that in a SUSY vacuum, (S h ) = 0. 22 However, that (S h ) 7^ can be 
shown explicitly by doing a one instanton calculation 

We thus face two choices, either {S h } ^ is a signal of the fact that 
the vacuum breaks SUSY or the relation S h = is modified by quantum 
corrections. There is overwhelming evidence that the second option is the 
right one. In particular the Witten index [32] for pure gauge theories is 
never zero implying that SUSY cannot be broken even by non-perturbative 
effects. 

The fact that the chiral ring relation S h = is modified to S h = const, 
is actually not that surprising. Remember that quantum corrections cannot 
bring new operators into existence but can modify the relations between 
existing ones by the introduction of the holomorphic scale A. The classical 
relation should be recovered by letting A — > 0. Simple dimensional analysis 
fixes S h oc A 3/l and the instanton computation reviewed in [23] (see also 
|24| ) allows also the determination of the proportionality constant C(G). 
We recognize in 3h the coefficient (3 of the one loop beta-function of pure 
SYM. From (|H4[) we can see that (S h ) is thus proportional to the exponential 
of minus the one instanton action 87r 2 /g 2 . Hence eq. (|145|) is changed into 

S h = C(G)A 3h in the quantum chiral ring. (148) 

So far we considered pure gauge theories. As it is going to be important 
later, let us consider an example of gauge theory with matter, say, a U(N) 
gauge theory with fields $ in the adjoint and Q, Q in the fundamental and 
anti-fundamental. Relation (|14()|) implies that, as far as the classical chiral 
ring is concerned, the fields W a and <3? commute inside various singlets and 
that W a Q = QW a = 0. Thus the generators of the chiral ring are of the 
following type jSSl El 01 : 

tr$ fc , trW a $ k , ti W a W a $ k , Q<f> k Q. (149) 

Classically, they are subjected to the relation analogous to (|145j) for S and to 
the usual relations stemming from the possibility of rewriting higher traces 
as functions of the lower ones. For instance, tr $ iV + 1 can be expressed as a 
polynomial in tr tr & N . 

When matter fields are present, there are also relations that use the dy- 
namics of the theory and depend on the particular choice of superpotential. 

22 By factorization one would also have (S) = 0. 
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Let $ r denote a generic field content, with r labeling the various irreducible 
representations of the gauge group. Let F r (W a , $>) be a combination of the 
fields that also transforms in the representation r. (Since we are working in 
the chiral ring, two such expressions that differ by a chirally exact operator 
can be identified.) The classical equations of motion for the matter fields 
are 

Contracting with F r to make a singlet this implies the classical relation 
0W($) 



-F r (W a ,<&) = in the classical chiral ring. (151) 



By now, it should not be surprising that also this classical relation is modified 
quantum mechanically. This is due to the so called Konishi anomaly which 
we now turn to discuss. 



6 The Konishi anomaly 

The Konishi anomaly is the SUSY extension of the usual chiral anomaly 
and it proves extremely useful in solving for the F-terms of the theory. 

In its original form it is based on the observation that (contrary to 
what we saw for chiral operators) the SUSY transformations of a generic 
composite operator are not always the ones that we would naively obtain 
applying the transformation to each component field. Consider again the 
example of SQCD, with its two chiral superfields Q and Q transforming in 
the fundamental and antifundamental representation of the gauge group. 
We saw in the previous section that the composite (chiral) operator qq con- 
structed from the lowest components obeys [Qa, qq] = as one would naively 
expect. 

If we drop the requirement that the composite operator be chiral, there 
are many other ways to build a gauge singlet even using just the components 
of one superfield instead of two. Of particular interest is the composite 
operator ip a q constructed only with the first and second components of Q. 
Naively, one would write: 

{Qd,#M = {Qa,r}q-r[Qa,q] 

dW 

= 2fq + = -2 — — q classical (152) 
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If l|152|) was true, the RHS could not get a v.e.v. in a SUSY vacuum. For 
instance, for W = mQQ one would deduce that (qq) = 0, contradicting 
many well established facts about SQCD. 

To obtain the right formula, we must once again regularize: 

^qix) = lim V^(x + e/2)(l - iefAJx) + ... )q(x - e/2). (153) 

Now in the SUSY variation of ip a q there will be an extra piece coming from 
the SUSY variation of A^: 

- -■ dW 1 

{Qd,#M = "2^g - -7= lim^ s + e2)ef i a? a X a (x)q(x - e/2) (154) 
oq V2 <^o ^ 

For our purposes it is enough to think of the gluino A as a background field, 
although this is not necessary. The extra piece seems to vanish as e — ► 
but now we have a singularity in the OPE coming from the presence of a 
term of the type y/2 J d^x q\ip in the Lagrangian. If one "brings down" such 
an interaction in the path integral one obtains the contribution (written for 
simplicity in the Abelian case) 

hme^^A Q (x)/'dV(0|r(^(x')^(x + e/2)) \0)\? (x')(0\T(q(x - e/2)q(x'))\0) 



x lime^^ a A a (x)e 4 ^^A /3 (x)^ oc X 2 (x). (155) 



More precisely, keeping track of all numerical factors and generalizing to 
the non- Abelian case one gets jS] 

{ Gd, ^q} = 2 (jq + tr X a \ a \ . (156) 
As usual, the same relation can be expressed in terms of superfields as 

This is so because we require the composite object Qe v Q to be a superfield, 
i.e. to contain a supermultiplet of fields properly transforming into each 
other by a SUSY transformation. For instance: 

qq S ^ q Q J q + _L tr (158) 
where the rightmost entry is thus the lowest component of D 2 Qe v Q. 
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In the application of the Konishi anomaly to the recent developments, 
it is sometimes necessary to consider a generalization of formulas (|156[) and 
(J157|) ,8! that we now turn to discuss. Consider chiral matter <3? in an 
arbitrary representation of the gauge group. Eq. I|157|) can be understood as 
the anomalous Schwinger-Dyson equation coming from varying $ — > $ + 77$ 
(77 <C 1) in the path integral. The classical piece comes from the variation 
of the classical action whereas the anomaly comes from the non-invariance 
of the functional measure. 

One can generalize the transformation to 

$->$ + 7/F($,W a ), (159) 

where F is a holomorphic function transforming in the same representation 
as <&. The classical piece is straightforward: 

D 2 [^e v F^=-^F classical. (160) 

As for the anomalous piece, the only additional complication is in keeping 
track of the color indices. Let r, s, . . . denote the indices of the representation 
in which $ transforms and let us write the gluino field as a matrix in that 
representation. The anomaly is, schematically 23 : 

lime^ r A*F s oc lime J i() p ) X s r X p q (F s ft), (161) 

where the spacetime dependence is as in Q155JI . Using 

BF - BF 

(Vrp P )<xS r p and (F.<ffl) oc -^(M q ) oc ^-Sf, (162) 

yields 

BF 

limeTp r \ s r F s oc A%— -. (163) 

q dct) q 

The generalized Konishi anomaly is thus: 

s 2 (*W)=-f^-^»^;fg. (164) 



23 There can also be corrections from chirally exact terms that would contribute to the D- 
terms, but the non-renormalization theorem prevents additional perturbative contribution 
to the F-terms. 
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The use of (|164|) is mainly in the context of the chiral ring, where we may set 
all chirally exact terms to zero, including the LHS of (|164|) . In this context 
(|TdH) reads: 



dW 



F r = - 



1 



W as r WJ q 



r dF s 



in the quantum chiral ring 



(165) 
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In particular, these relations hold for v.e.v.s in a SUSY vacuum, and they 
are often enough to solve for the F-terms of the effective theory. Specifically, 
they can be used to prove the Dijkgraaf - Vafa conjecture as we now discuss. 

7 Glueball superpotential via the Konishi anomaly 

In this section, we will exploit a relationship which exists between the ef- 
fective superpotential and the SUSY vacuum expectation values of chiral 
operators to determine completely the former. The relation was discussed 
in Section 3 and reads: 



This relation was derived when we realized that integrating out and back in 
degrees of freedom was technically a Legendre transform. Alternatively, one 
can derive the above expression directly by expressing W e g in terms of the 
path integral over the elementary fields, weighted by the action including the 
tree level superpotential Wtree = X) ^k^k ■ Then considering the coupling 
Afc as a source term for the operator the relation (|166j) follows. 

The usefulness of eq. (|166[) is that if we are able to determine (Xj.) as a 
function of the couplings A^ and of S we can then obtain W e g up to a function 
independent of the coupling. The Konishi anomaly and its generalizations 
provide us with just a way to determine such expectation values. 24 

Let us see how this works in a simple example, where the non-generalized 
Konishi anomaly is sufficient to solve for a non-trivial W e g |36| . We take 
SQCD with gauge group U(N) and, for simplicity, only one flavor 25 . 

The only gauge invariant we can build with this matter content is a 
single meson superfield, M = Q a Q a - A tree level superpotential that one 
can write is the following: 



That the Konishi anomaly plays an important role in the Dijkgraaf - Vafa conjecture 
was first emphasized in I35|. 

25 The C/(l) factor is unimportant as it decouples in the IR. 




(166) 



Wtree 



mM + AM 2 . 



(167) 
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The quarks are thus massive, and the second term introduces a non-trivial 
interaction. The latter is non-renormalizable, but it is straightforward to see 
that (|167|) can be obtained by writing the renormalizable superpotential: 

W tree = mQQ + |tr $ 2 + gQ<S>Q, (168) 

and then integrating out the adjoint field <3?, assuming that \/j>\ 3> \m\ and 

2 

identifying A = — 

Let us now write down the Konishi anomaly relation for, say, the U(1)q 
rotating the superfield Q. This is basically the computation leading to (|157|) . 
Using ()167l) . the relation in a SUSY vacuum becomes: 

m(M) + 2A(M 2 ) = S. (169) 

By using factorization, we can reexpress (M 2 ) = (M) 2 and then solve (|169|) 
for (M): 

The two signs in (jl70[) have a clear meaning in the classical limit, where 
S — * 0. For the positive sign, the v.e.v. tends to zero, and we are thus in 
a situation where classically the U(N) gauge group in unbroken. For the 
negative sign, the v.e.v. tends to a finite value, and thus U(N) is broken to 
U (N — 1). Note that when the coupling A is sent to zero this second classical 
vacuum is sent to infinity and disappears. 

We thus have two partial differential equations for W e g: 

^ = <M(S,m,A)>, ^ = (M(S,m,X)) 2 - (171) 
Using the expression (|17U|) . we obtain the solution: 



w * ~£ ^f^W s+ s^ + sv*^ f^sj + CIS), 

(172) 

where of course the coupling independent piece of W e g is still undetermined, 
and we have used the holomorphic scale A to compensate for the dimension 
of m. This can be done because a rescaling in A can always be absorbed in 
C(S). 

The coupling independent piece C(S) can be determined by using low- 
energy information on the theory. Consider simply the positive branch of 
(|17()|) . and take the coupling A to zero. Here we should recover the low energy 
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physics of SQCD with a massive flavor. We know what this is, namely 
pure SU (N) SYM, and we also know it is described by the Veneziano- 
Yankielowicz superpotential (|75|) . By thus imposing that (J172p with the 
positive signs reduces to (|75|) when A — > (for this, we have to match the 
holomorphic scale A of the pure SYM theory to the scale A of the theory 
with one flavor, A 3N = mA 3W_1 ) we obtain: 

W cS = Wl-logJ^+Slog^-l (173) 

m 2 , m 2 r 8A~ F 8A~\ 

r±— r\/H -S + Slog -±-\ l-\ =S . 

8A 8A V m 2 y2 2 V m 2 j 

We have thus determined the full effective superpotential by using the Kon- 
ishi anomaly (|169|) to determine the v.e.v.s, which have allowed us to solve 
for the coupling dependence of W e g. The remaining piece is then deter- 
mined by taking specific values of the coupling constants which make the 
low-energy physics more familiar. 

It is worth noting that the same result (|173|) can be obtained through 
a matrix- vector model based on the tree level superpotential (|168|) |37j . 
There one can show that only the graphs with the topology of a disk give a 
non-trivial contribution, the boundary of the disk being composed of quark 
propagators. The result of this perturbative computation is a power series 
in S, which sums to the exact expression (|173j) . In this latter approach the 
VY piece of W e g has to be added to the result. 

We now turn to the example discussed in |Sj, where the generalized 
Konishi anomaly relations will be needed. We consider a U(N) gauge theory 
with a matter field in the adjoint representation and a generic tree level 
superpotential: 

n 

^t r ee = Etfr tr$fc+1 - (174) 

k=0 

The classical vacuum structure is determined by the extrema of W tree (z), 
seen as a holomorphic function of a complex variable. If we assume that the 
n extrema are all distinct: 

n 

w Le( z ) = 9n Y[(z - <h), di ^ cij for i ^ j, (175) 

i=\ 

then in the classical vacuum the eigenvalues of $ will partition into groups 
corresponding to every extremum and the fluctuations around these vacua 
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will all be massive. If we have Ni eigenvalues equal to aj, the classical gauge 
symmetry is spontaneously broken in the following way: 

n 

U(N) - HU(Ni). (176) 

i=i 

At the quantum level, we expect for every U(Ni) factor that at low enough 
energies the SU(Ni) piece will confine, while the £7(1) piece will be in a free 
Coulomb phase. For every factor, we thus define the effective fields: 

Si = -^trW?W ai , (177) 
wf = ^trWf. (178) 

The quantity that we want to determine is W e s as a function of the couplings 
gk and the effective fields Si and wf. Notice that the trace in 1)1 77|) is over 
the full U(Ni) subgroup and not SU(Ni) for later convenience. 

We now make a crucial remark that simplifies the discussion. The overall 
£7(1) of the full U(N) gauge group is free, since the matter field is in the 
adjoint representation. This also means that there are no Yukawa couplings 
involving the gaugino relative to this £7(1), and thus its action is also one of 
a free fermion. We thus have the possibility of shifting the £7(1) gaugino by a 
constant Weyl spinor without changing anything in the theory. In superfield 
terms, this symmetry of the theory is implemented as: 

W a -» W a - 4TT Xa lNxN. (179) 

When the gauge symmetry is broken as in ()176|) . the symmetry still acts as 
an overall simultaneous shift of all the gaugino superfields, W a i — > W a i — 
A-KXa^NixNi, so that the effective fields transform as: 

1 

Si -> Si + x a w ai - -NiX a Xa, 

wf -> wf-N iX a - (180) 

The effective superpotential W e s(Si,wf ,gi~) has to be invariant under the 
above transformation, since it cannot depend on the decoupled overall £7(1). 

A clever way to constrain the form of W e g using the invariance under 
(|18U|) is to apply to this fermionic symmetry the same tricks used with su- 
persymmetry in Section 2. Namely, consider the invariance as a translation 
in an auxiliary Grassmann variable vp a . Then build a "superfield" using 
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ip a such that the translation tp a — > ip a + x a implements ()18U|) . Finally, an 
invariant function is trivially built by writing the integral over ■i/'-superspace 
of any function of the ^-superfield. Hence, we define: 

A manifestly invariant effective superpotential is then given by: 

W DV = - Jd 2 ^(S h g k ), (182) 



which, after performing the Grassmann integral, becomes: 

pf-p 1 pp, 



To compute the derivatives we can consider T as a function of the Si , that 
is the lowest components of the Si. 

Note that ()183|) is quadratic in wf, a result which could also be obtained 
in the perturbative set up of Section 4 by observing that only planar di- 
agrams contribute. We should stress here that this feature appears only 
when Si is the glueball field obtained by tracing over U(Ni) and not over 
its SU(Ni) subgroup. If we were to consider the latter, then W e ff would 
contain higher powers of wf. 

Let us anticipate now that the main result of this section will be to 
identify the function J-(Si,gf.) with the planar free energy of the related 
matrix model, as eq. (|124jl at the end of Section 4 suggests. To see how this 
comes about, let us derive equations for this function T . 

First of all, the coupling dependence of the effective superpotential is 
fixed by the equations: 

*Wk-/^tr*H-A ( 184 ) 



dg k \k + l 
Using eq. p82j> . the above equation reads: 



/ 



rf2 ^ (s " ffi) = -(iTl"' +1 >' <185) 
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We can write an equation directly for the integrand T if we write the RHS as 
the highest component of a t/'-superfield, generalizing ()181j) . We thus have: 

For convenience, the equations can be written setting tp a = without loss of 
information, since from the lowest components the higher ones are generated 
acting as in (|179|) . The equation we will thus take as defining the coupling 
dependent part of the function T is: 

The next step is to introduce a generating function for the operators on the 
RHS of the above equation. It reads: 

*« = -^( tr ^) S -^t^i^W^). (188) 



32tt 2 \ z - $ / 32?r 2 ^ z k 

Note that in the above we are not really paying attention to the ordering 
of the W a with respect to the <£, the reason being the relation in the chiral 
ring \W a ,$\ = 0. 

The expression (|188|) also gives us a more precise, gauge invariant, def- 
inition of the gaugino condensates Si. Indeed, suppose that the analytic 
structure of the function R(z) is such that it has n cuts over the z-plane 
(we will shortly see how this arises) . The cuts can be seen as a sort of quan- 
tum resolution of simple poles. By defining the contour as the one going 
around the i-th. cut, we define: 

Si = — <p dz R(z). (189) 



2tti jq. 

Classically, the above relation is understood as follows. If the contour d 
encircles the eigenvalue aj, then for any matrix M the integral projects on 
the eigenspace corresponding to this eigenvalue. Indeed, having diagonalized 
we have: 



i f , + m i r A M mm 

az tr = dz > 

2tti J Ci z-Q 2m J Ci ^ z - 



= 2^ 



""•>Vi m= i -rrn 

M mm = tvP t M = tvMi. (190) 
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Note that we also have trPj = JVj. Since projectors should not receive 
quantum corrections, the result of the manipulation above carries over to 
the quantum level. 

We are now after a way to solve for the generating function (also called 
resolvent, in reference to the matrix model analogy) R{z). We are going to 
see that it is precisely the generalized Konishi anomaly which provides us 
with a closed equation for it. 

Let us first of all specialize eq. (|165j) to the case of a matter field in 
the adjoint. The LHS simply becomes a trace. As for the RHS, we have 
to make explicit the adjoint action of the gaugino superfield by writing 
W ar s = iW ap f pr s , with all indices in the adjoint representation. Then the 
two structure constants are traded for two commutators, and the generalized 
Konishi anomaly reads: 



<tr = £ (\w 



w. dF 



(191) 



We now specialize to a particular variation of <£, namely we take: 



1 W a W n 
62tt z z — <P 

To see what the RHS of ()191|) looks like with this variation, we develop in a 
power series as in (j!88[) : 



d (W a W n \ k ^ 1 d 



\n\k 



d<£>\ \ z-$ J, ^ z n+1 <9$ 

J 1 n=0 



Z-/ z n+l ^ 
n=0 3 

00 n— 1 

E E ^ TT (w a w a ^)H^ p - 1 )i 

n=l p=0 

00 00 

E E -^ ri (w a w a $nH$ n - p - l )i 

p=0 n=p+l 
00 00 ^ 

E E i^r^(^^* p )* (* m )?' 



p=0 m=0 



$ / , \ Z - $ / ; 



(193) 



Implementing now the chiral ring relations, which forbid traces with more 
than two W a , we see that only one term of the (anti)commutators in (|191|) 
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is non-vanishing, namely when the four W a distribute themselves by two in 
each trace. We thus get the following basic relation in a SUSY vacuum: 

-=.H-^))-w((^)}' (I94) 

At this point, we can crucially use the factorization property of the correla- 
tion functions of chiral operators to write the RHS as R(z) 2 . Concerning the 
LHS, we can add and subtract W( ree (z) to W{ Tee (^) so that we can write: 

"3^ ( tr (<-(*)ttf) ) = + ( 195 ) 

where we have defined: 

fn-i(z) = ± (tr {{WU(z) - < ee ($))^^ ^ . (196) 

It is easy to convince oneself that f n —i(z) is a polynomial of degree n — 1: 
fn-i(z) = Y^k=o fk zk ■ Indeed, by expanding W^eeC 2 ) m the numerator 
around any eigenvalue of <£, one can see that f n -i{z) has no singularities. 
Moreover, we see that at large z, the LHS of ()195j) goes like 1/z. Thus 
f n -i(z) must cancel all the non-negative powers in the product W( ree (z) R(z) , 
which go at most as z n ~ l . At this point the n coefficients of f n —i(z) are 
unknown complex parameters, to be related to more physically significant 
quantities shortly. 

We have thus finally obtained the equation for R(z). It reads: 

R(z? = Wl ree {z)R{z) + \f n ^(z). (197) 
Solving the above quadratic equation yields 

= \ (<ee(*) " V^treeW' + Zn-lW) , (198) 

where the sign has been chosen in order to get the right 1/z behavior at 
infinity. The above solution gives us an expression of R{z) (and thus, ex- 
panding, of all the relevant v.e.v.s) in terms of the couplings gj. and the 
coefficients fk of the polynomial. 

In order to have a useful solution for R(z), namely one which is a function 
of <7fc and the Si, we only have to use the relations (|189|) to get expressions 
of the Si in terms of the fk and g^. By inverting these relations, we finally 
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obtain an expression R{z) in terms of gk and S{ that can be used to integrate 
for !F{Si,gk). Fortunately, in this case we do not need to perform this 
computation explicitly because the implicit expression (|189|) is all we need 
to show the correspondence with the matrix model. 

We now pause to mention what should be done in order to perform a 
more complete treatment of this theory. One defines two more generating 
functions: 

-w = i( tr r^)' rw = ( tr ^}' <199) 

which can be seen as higher components of R{z) in a ^-superfield. Anomaly 
equations for them can also be written, which can be solved linearly. New 
polynomials with undetermined coefficients like f n -\(z) arise, but again 
these coefficients can be solved for the physical quantities: 

In principle, it is the expansion of T(z) that is needed to obtain W e ff, but 
as we have seen, the ip translational symmetry makes it possible to consider 
only the resolvent R(z) to fully solve the problem. This is no longer true 
for some theories with a different gauge group and/or matter content, for 
which one cannot write the fermionic if) symmetry. One can still write 
generalized Konishi anomaly relations that can be solved, but in this case 
one has eventually to get the expression of the analog of T{z) in order to 
obtain W e s- 

What we are going to show in the following is that there is an even more 
direct way to compute J-{Si,gk) through a bosonic matrix model computa- 
tion. To prove this, we are going to show that exactly the same equations 
for R{z) and T arise within the matrix model. 

Let us define the free energy T of the matrix model whose potential is 
given by W tr ee- 

e~W f = JdM e"T "^W, (201) 

where M are hermitian N x N matrices and we take N to be large. Note 
that N is not related to the N of the gauge theory. For the moment we 
may also think S as unrelated to S although a relation between these last 
quantities will arise in the following. 
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First of all we derive differential equations for T by taking derivatives 
with respect to the couplings on both sides: 



dT Si i k+1 



dgk N \ k + 1 



tr M + ) . (202) 



In the large N limit, both sides of the above equation will have a genus 
expansion, with the planar contribution being the leading one. (With this 
normalization they both tend to a finite limit.) 

As in (|188|) . we introduce a generating function for the v.e.v.s on the 
RHS of (|2()2j) . This is called the matrix model resolvent: 

We are now left with the task of deriving an equation for the resolvent. 
Similarly to the Konishi anomaly, we can obtain such an equation from 
the Ward identities implied by reparameterization invariance of the matrix 
model free energy. Namely, the integral in (|2()lj) should be invariant under 
a simultaneous reparameterization of the matrix M — > M + 5M both in the 
measure and in the potential. Let us take the following variation: 

5M = 7] — — , T) < 1. (204) 



The measure changes as: 

/ imp. 



( d5M j \ 
dM dM 1 + %- . (205) 



dM? 

Now, using the same trick as in (|193|) . we have that: 



9 ( 1 V ^-^1 • (206) 



8Ml \z-Mj l V z-M 

On the other hand, the variation (|2U4[) acts on the potential by bringing 
down its derivative, so that the Ward identity reads: 

tr EkH>\. (207 ) 

z-M I v ; 
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This expression can be turned into an equation for R(z) again using factor- 
ization, but this time the reason we can factorize is very different from the 
one in the gauge theory context. Indeed, we can write: 

only in the planar limit, that is, keeping only the leading order in N for 
both expressions. 

This can be understood as follows. Consider the correlation of two op- 
erators (l/iV)trM 4 as depicted in Fig. [Tj The diagrams contributing to 




Figure 7: Factorization - the disconnected two-point function of the operator 
(l/iV)tr M 4 is finite while the connected one is suppressed by a factor 1/N 2 . 

the correlator have the same form of vacuum diagrams, only this time at 
the point where the operator is inserted we have a factor (1/iV) instead of 
N as there would be for an ordinary vertex. Thus we see that connected 
diagrams with one insertion are finite while those with two insertions are 
suppressed by (1/iV 2 ) leading to the abovementioned factorization property. 
The mechanism leading to factorization here is thus totally different from 
the gauge theory context, where supersymmetry and the properties of chiral 
operators are crucially used. 

As for the RHS of (|207|) . we apply the same addition and subtraction 
of W( rcc (z) to write the final equation for the resolvent, known as the loop 
equation, as: 

R(z) 2 = W( tee (z)R(z) + \fn-i(z), (209) 

with: 

also an unknown polynomial of degree n — 1. We thus see that the loop 
equation for R(z) is the same as eq. (|197j) for R{z) and thus that R(z) 



58 



has the same analytic structure as R(z) in the gauge theory, namely it will 
display a total of n cuts. 

As before, the coefficients of the function f n -i( z ) can be reexpressed in 
terms of the following variables: 

Si = I dz R(z), (211) 

where the contour Ci goes around the z-th cut 26 . Using the definition ()2U3|) 
of the resolvent, and applying the same reasoning as in l|190|l . we obtain: 

Si = — I dz ^ I tr — !— \ = -2-JVi. (212) 
27ri Jq. N \ z — M I N K ' 

The Si are proportional to the fraction of eigenvalues populating the i-th 
cut 27 . 

We are now finally ready to state the correspondence between the matrix 
model and the gauge theory computation. We already have the same poten- 
tial Wtree(z) entering the equation for the resolvent. We now identify also 
the two polynomials of degree n — 1, that is we impose f n ~i(z) = f n -i(z) 
so that the coefficients entering the loop equations are all the same. We can 
then trivially identify the two resolvents: R(z) = R{z). This in turn implies 
Si = Si. Lastly, expanding the resolvents we also obtain identical equations 
for the free energy on one side and for the function T on the other, so that 
we identify T{Si,g k ) = T(S i: g k ). 

This is thus the core of the correspondence between the computation of 
the effective superpotential on the gauge theory side, and the free energy 
of the matrix model. Using a mixed notation, and recalling the relation 
between T and \Vdv, we have that: 

a-ir 

Wdv = ^ N ^ ( 213 ) 

i 

where we concentrate on the piece of Wdv relevant purely to the glueball 
super fields. 

There are some important remarks one has to make on (|2L-ij) . The first 
is to state once again that the iVj are not related to the iVj, the former being 



Note that in the hermitian matrix model context, the cuts are on the real axis and 
thus the Si are real, unlike the gauge theory context where both can be generally complex. 
The relation between the quantities on the two sides will be done by analytic continuation. 
27 This last relation has no analogue in the gauge theory. 
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finite and the latter being taken to infinity in order to extract the planar 
contribution to the free energy (and to be able to use the factorization in 
lHU)). 

The second remark is that once the correspondence is derived as before, 
one can turn to any other alternative derivation of T on the matrix model 
side in order to compute it, as for instance presented in Appendix B. Then, 
when the result is stated in terms of the filling fractions S% and the couplings 
gk, the gauge theory effective superpotential is obtained using (j213j) . 

The last comment on (|213|) is that one could take it at face value and 
apply it even to any part of the matrix model free energy T which depends 
only on the filling fractions Si and not on the couplings gt- That would 
imply the possibility of deriving on the gauge theory side the effective su- 
perpotential also for the pure gauge low-energy dynamics. Keeping in mind 
that all the arguments of this section do not apply to this case, one can still 
perform a simple computation and find, very similarly to (jl26[) . that the 
Veneziano-Yankielowicz superpotential is reproduced. This fact still needs 
a better understanding in gauge theoretic terms. 
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A Superspace notation 

In this appendix we collect the notation and conventions used throughout 
this paper. 

The Minkowski metric is rf- v = diag(+l, —1,-1,-1). 
The gamma matrices 7^ satisfy the Clifford algebra {7^,7^} = 2//^, 
and we always use the Weyl representation 

f=(L f), (214) 



where 



at* 



° ~ 1 -1) ' a = ~ a = [1 oj' 



and the indices are read, row by column, as (tr^)^ and {a^) aa . For instance, 
(^ 2 )i2 = 

The generators of the Lorentz transformations are 

= \(o^ a a v h -a v « a o»\. (217) 

Dotted and undotted indices describe two component Weyl spinors of 
opposite chirality and we raise them and lower them using the totally anti- 
symmetric tensors £21 = e 12 = 1 and = e 12 = 1: 

V> a = e°%, ^ a = e a ^, V^ = e A %, $& = e & ^- (218) 

Complex conjugation changes chirality: (ijj a )* = vpa and when two pairs 
of indices are not explicitly written, the following summation conventions 
are understood: 

= fx« and = $aX a - (219) 

So far everything is as in except for the signature of the metric. We 
now define the square of a spinor as 

^ 2 = ^ Q Va and $ 2 = ■ (220) 
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(This is done as to eliminate many of the powers of two appearing explicitly 
otherwise.) 

Superspace has coordinates x^,9 a 1 a and the fermionic derivatives are 
defined as 

8 a = — and 8& = --t, (221) 
d6 a d9 a v ' 

so that d a 0P = bl and d & $ = tfg. 

The supercharges and superspace derivatives are defined as follows: 

Qa = d a -~a» & e% and Q« = 4-VC^ ( 222 ) 

D a = 8 a + and D a = d a + ^9 a a^. 

The only two non-vanishing anticommutators are 

{Qa,Qp} = —{D a , Dp} = -ia^d, = -a^ $ P,. (223) 

Indices on the derivatives D a and Da are raised and lowered using the 
same conventions as in (|218|) but we define the square of the fermionic deriva- 
tives with an extra minus sign as compared to (|220f> : 

D 2 = -^D a D a and D 2 = -^D a D 6: , (224) 

so that: D 2 9 2 = D 2 9 2 = +1. The fermionic integrals are just another 
notation for the derivative: J d 2 6 6 2 = J d 2 9 9 2 = 1. It is straightforward 
to check that these integrals are invariant under the translations 9^9 + ^ 
and # ^ + 

All the identities involving the a matrices described in the Appendices 
A and B of Wess and Bagger 1.1 j are still valid by letting r]^ — ► —r]^ u and 
recalling that 99 = 29 2 and 99 = 29 2 in our notation. 



B One cut solution to the cubic matrix model 

In this appendix we show how to compute the leading (planar) contribution 
to the free energy T of the one cut cubic matrix model jlUj , thus completing 
the computation of Section 4. To be slightly more general, we present the 
computation for a generic potential W and only at the end we substitute 
W(z) = fz 2 + f z 3 . 
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We consider the following matrix integral for a generic potential W(M) 
of a random hermitian matrix M (here we write the trace explicitly) 

e -S^ = |d%e-f trW ( M ', (225) 

where T is the planar free energy, the dots denote non-planar (sub-leading) 
contributions and with respect to Section 4 we have already substituted for 
e = S/N. 

Since M is hermitian, there exists a unitary matrix 0, such that 

M = Q^AU , (226) 

where A = diag (\\, ...,Aa) is a diagonal matrix with the real eigenvalues 

Aj of M as entries. The measure in the new variables becomes d N2 M = 
d N Ad N2 ~ N 0, A(A) 2 with A(A) 2 being the Jacobian of the transformation. 

This Jacobian can be determined in the following way. The metric on 
the matrix space can be written as 

ds 2 = tr dMdM 

= tr (anUn + rtdAn + ftUdft) 2 

= tr(dA+ [A,d^ f ]) 2 
= tr(dA 2 + [A,df^ f ] 2 ) 

= dA 2 + ... + dA 2 ^ + ^(A i -A j ) 2 |d^y 2 . (227) 

The Jacobian is just the square root of the determinant of this metric 

A(A) = Vdet G = JJ(Ai - Aj) 2 . (228) 

i<j 



This result can be inserted into eq. (|225|) 

jd^Me-i trW(M) = Jd**-*nd*\ n(A, - \ 3 fe~^ W ^\ (229) 



The integral over the unitary matrices gives an overall constant that we 
absorb into the measure leaving 

/d* 2 Me-f W(M) = jdhe-^ W{Xk)+ ^ lo ^- X ^. (230) 
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Solving this integral is impossible for finite N. Instead, we use the steepest 
descent method for the evaluation, which corresponds to taking the large N 
limit, thus singling out the planar diagrams. 

To leading order in N, the free energy is given by 

-f^ = -f E w ^ + E lo § - U ( 231 ) 

k i^j 

where A& are the N solutions to the extremization problem given by the 
simultaneous solution of the following N equations, labeled by k = 1, ... N. 

= -^ W ' (Ajfc) + 2^_L-. (232) 
S k+o k j 

To solve this set of equations for large N one goes to a continuous de- 
scription of the eigenvalues by defining 

4 -» x G [0, 1] 
N 

V -► N [ dx (233) 
Aj ->■ A(x) 

With these identifications the free energy and the extremum equation be- 
come 

T = S [ dxW{X(x)) - S 2 [ dx [ dylog\X(x) - X(y)\ (234) 
Jo Jo Jo 

= -W^^sfa^L—, (235) 

where the integrals are to be thought of as principal values. Note that all 
dependence on N has disappeared and that we traded a set of TV algebraic 
equations for one integral equation. 

To proceed, introduce the density of eigenvalues p(X) as 

Sdx = dX p(A) (236) 

which is non-negative, has support within an interval [a, b] (to be deter- 
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mined) and is normalized such that 



f dAp(A) = S. (237) 

J a 

In terms of the eigenvalue density eqs. ()234j) and (|235j) read 

T = ! d\p{X)W{\)- f dA' f dAp(A)p(A')log|A-A'| (238) 

J a J a J a 

= -W'{\) + 2 f dA' S^Ij . (239) 

J a A — A 

Integrating (|239|) in the two intervals [a, A] and [A, b] and combining the 
results yields 

f dA' p(A') log |A - A'| = j (2 W{\) - W(b) - W(aj) 
+i J dA'p(A') log(6 - A') (A' - a) (240) 

Eq. ()238|) then becomes, using the normalization of the eigenvalue density 
(12371) 

•T 7 = \j\\ P {\)(w{\)-Slog{b-\){\-a)) 

+ ffw(6) + W(o)). (241) 



4 

To find the eigenvalue density define the resolvent 



00 



/?(:)-/ <1A^, (242) 



00 



which is an analytic function of z. The singularity of the resolvent is a cut 
along the real interval [a, b\. Across the cut the resolvent has the following 
behavior 

R{\ + ie) + R(X -ie) = 2 f dA' = W'{\) (243) 

J a A - A 

i?(A + ie) - R(X - ie) = jd\' = -2vrip(A) . (244) 



28 We focus here on the so called one cut solution, i.e. the one with all eigenvalues 
grouped around an extremum, corresponding to an unbroken gauge group. 
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The integral in Q243JI is a principal value and the one in (|244j) is a residue 
as can be easily seen by drawing the contours. The idea now is to use ()243)1 
to find the resolvent and then 1)244(1 to extract the eigenvalue density p. 

One particular solution to eq. (|243|) is, of course, the regular solution 
Rreg( z ) = ^W(z), but this cannot be the full story because it does not have 
the right asymptotic behavior at infinity. To R reg (z) we must add a generic 
solution to the homogeneous equation 

R ain g(\ + ie) + R ain g(\-ie) = 0, (245) 

necessarily singular along the cut and chosen so that the full solution has 
the correct large z behavior following from 1)242)1 : 

R(z) = R sing {z) + R reg {z) ~ | + O (J^J > z - . (246) 

Across the cut we need a sign change which is implemented by 

Rsing(z) = P{z)^{z-a)(z-b), (247) 

where P{z) is a regular polynomial with degree two less than the degree of 
W(^) in order to cancel the various powers of z. The coefficients in P(z) 
are given implicitly by the condition (|246|) . 

The sought eigenvalue density p(X) is then, using eq. 1)244)1 

p(X) = —R smg (X), (248) 

TTl 

where the sign of the square root is chosen to give p non negative. Having 
determined the eigenvalue density, the free energy can be calculated by 
evaluating (J2HJ). 

Up to this stage the potential W(M) is generic. To illustrate the tech- 
nique and to give the full solution to the problem of Section 4, we now 
specialize to a cubic matrix model. Notice that although the integral with a 
cubic potential is not convergent, each term in the perturbative expansion 
is, and this is all we need. Let us thus take, as in Section 4 

W{M) = —M 2 + -M 3 . (249) 
2 3 

Requiring that there are no 0(z 2 ) and O(z) terms in R(z) in the large 
z limit, we get: 

1 

P(z) = — (2m + g(a + b) + 2gz) . (250) 
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Furthermore, by canceling the constant term O(z ) and requiring the pole 
to have residue S, we obtain conditions on the end points a and b of the 
interval where the eigenvalue density has support 

= 2g(a + b) 2 + g(a - bf + 4m(a + b) (251) 
165 = {a - bf {g {a + b) + m) . (252) 

Thus, from eq. 

1 



p(\) = — (2g\ + 2m + g (a + 6)) a/(6 - A)(A - a). (253) 

47T 

The free energy is computed by using the eigenvalue density and evalu- 
ating eq. flMJ) 



- /dA [(2 5 A + 2m + g(a + 6)) y/(b - A)(A 

'71" J a L 



^X 2 + ^X 3 -Slog\(b-X)(X-a)\ 



+f(f(a 2 + 6 2 ) + f(« 3 + 6 3 )) (254) 



4 V 2 v ' 3 

It is convenient to introduce a parameter 



a = — (a + b). (255) 
2m 

Then the constraints (|251j) and l)252|) imply 

2-^rS + a{l + <t)(1 + 2tr) = 0, (256) 

where we chose the branch a = for g = 0. Evaluating the integral Q254|) 
keeping only terms that depend on g (through a) yields 29 : 

g,r(2 + tt, + 3^) + S2 

3 (1+ct)(1 + 2ct) 2 2 sv ; v ; 

The constraint (|256|) can be solved in a perturbation expansion in the cou- 
pling constant g and substituting into T we get: 

jr- 2 9 2 & 8 g 4 ^ 4 56g 6 S 5 512g 8 ff 6 9152g 10 5 7 | g 12 
3m 3 3m 6 3m 9 3m 12 5m 15 

The superpotential for the gauge theory is then given by the Dijkgraaf- 
Vafa formula (|124f) of Section 4. 

29 A constant term arising at g = can always be absorbed into the measure as we have 
done many times before. 
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C Bibliographical note 



As mentioned in the Introduction, the work of Dijkgraaf and Vafa has 
generated many other interesting developments that unfortunately did not 
find their place in our review. In this Appendix we make an attempt to 
summarize them hoping that this will help the reader in finding his /her own 
way through the literature. We apologize for any omission or imprecision 
that might occur. 

The most important omission is, of course, a discussion of the geomet- 
rical aspects and implications of the work of Dijkgraaf and Vafa. This is 
particularly embarrassing because this was the original route through which 
the equivalence with the matrix model was originally discovered! This sub- 
ject deserves its own review written by someone more qualified than us to 
do it justice - here we only refer the interested reader to some of the related 
literature [SHI that followed p. 

Within the realm of gauge theories, probably the most important omis- 
sion is the discussion of theories with extended SUSY. The work of has 
a natural extension to theories with N = 2 and N = 4 SUSY. In the N = 2 
context [HHl many of the results of Seiberg-Witten theory [H] have been repro- 
duced using matrix models. Various extensions to gauge groups other than 
U(N) and the addition of flavor hypermultiplets have also been discussed. In 
the A/" = 4 case 0U], much attention has been given to extracting S-duality 
and to the various deformations allowed by the presence of three chiral mul- 
tiplets in the adjoint representation. In particular, one can study either 
the Leigh-Strassler deformation [JJ or the relevant deformations caused by 
adding mass terms to some of the chiral multiplets leading to the so called 
N = 1* or N = 2* theories @2j. 

The new techniques to study SUSY gauge theories immediately call for 
their application to the study of various non-perturbative aspects of their 
dynamics. We greatly regret not having discussed the new results emerging 
from these investigations such as a much better understanding of the pa- 
rameter space of these theories [13] , (see also [H] ) , the analysis [IHj of the 
Argyres-Douglas points jlH] and that 07j of Seiberg duality [iH] . 

Another notable omission is the discussion of gravitational corrections 
|49j . In the presence of a non-trivial background for the gravitino, the 
analysis that led to the vanishing of non-planar diagrams is no longer valid 
because the gravitino can soak up additional zero modes. The non-planar 
contributions have also been checked to agree with previously known results 
with both a diagrammatic analysis and a computation based on the anomaly. 
This has also led to an interesting deformation of the chiral ring where 
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"Grassmann" operators are no longer nilpotent even at the classical level. 

Returning to gauge theory, the matrix model has been applied to the 
study of theories based on other gauge groups such as SO(N) and Sp(N) 
without |5L)| or with [SJ flavors as well as U(N) with symmetric or anti- 
symmetric representations. (The case of U(N) with "flavor" matter in the 
fundamental has already briefly been discussed in the main text jSZ], see 
also |52j.) The case of SU(N) with a number of flavors equal to the number 
of colors is particularly interesting because it allows for baryonic deforma- 
tions [5Hj . Various aspects of this model have been investigated in [54| . 

Moving on to more "exotic" cases, orbifolds and quivers yielding semi- 
simple groups have been studied in whereas multi-trace deformations 
have been discussed in [55] . 

Regrettably, the study of chiral gauge theories using these new techniques 
has not received much attention. These theories are interesting because 
they can [HSj and they do j^Z] break SUSY dynamically. Some examples are 
discussed in [3BJ, see also [55] . 

Still within the realm of applications to four dimensional gauge theories, 
a justification of the linearity principle of 2Q] from the matrix model has 
been proposed in jS^] and an attempt at deriving the VY term from a 
diagrammatic expansion has been made in jSUj. A thorough study of some 
of the more formal aspects of the matrix model has been done in the series 
of papers [61] . (See |H2] for a discussion of the gauged model and the role 
played by the ghosts and |k5] for some studies on supersymmetric matrix 
models). Finally, although the original results are strictly four dimensional 
as should be clear from the discussion of Section 4, extensions of some results 
to both lower [12] and higher [HE] dimensions have also appeared. 



References 

[1] R. Dijkgraaf and C. Vafa, |arXiv:hep- th/0208048 , 

[2] D. Amati, K. Konishi, Y. Meurice, G. C. Rossi and G. Veneziano, Phys. 
Rept. 162 (1988) 169. 

[3] K. A. Intriligator and N. Seiberg, Nucl. Phys. Proc. Suppl. 45BC 
(1996) 1 |arXiv:hep-th/9509066| . 

M. A. Shifman, P rog. Part. Nucl. Phys. 39 (1997) 1 
|arXiv:hep-th/9704lT4] . 



D. S. Berman and E. Rabinovici, arXiv:hep-th/0210044 



69 



[4] R. Gopakumar and C. Vafa, Adv. Theor. Math. Phys. 3 (1999) 1415 
|arXiv:hep-th/981113l"l . 

C. Vafa, J. Math. Phys. 42 (2001) 2798 |arXiv:hep-th/0008142| . 

F. Cachazo, K. A. Intriligator and C. Vafa, Nncl. Phys. B 603 (2001) 
3 |arXiv:hep-th/0103067| . 

R. Dijkgraaf and C . Vafa, Nucl. Phys. B 644 (2002) 3 
|arXiv:hep-th /0206255] . 

R. Dijkgraaf and C . Vafa, Nucl. Phys. B 644 (2002) 21 
|arXiv:hep-th/0207i06] . 

[5] F. Ferrari, Nucl. Phys. B 648 (2003) 161 |arXiv:hep-th/0210135| . 

[6] N. Seiberg and E. Witten, Nucl. Phys. B 426 (1994) 19 [Erratum-ibid. 
B 430 (1994) 485] arXiv:hep-th/9407087| . 

N. Seiberg and E . Witten, Nucl. Phys. B 431 (1994) 484 
|arXiv:hep-th/9408'099] . 

[7] R. Dijkgraaf, M. T. Grisaru, C. S. Lam, C. Vafa and D. Zanon, 
|arXiv:hep-th/02110T7l 

[8] F. Cachazo, M. R. Douglas, N. Seiberg and E. Witten, JHEP 0212 
(2002) 071 arXiv:hep-th/0211170| . 

[9] K. Konishi, Phys. Lett. B 135 (1984) 439. 

K. Konishi and K. Shizuya, Nuovo Cim. A 90 (1985) 111. 

[10] E. Brezin, C. Itzykson, G. Parisi and J. B. Zuber, Commun. Math. 
Phys. 59 (1978) 35. 

[11] J. Wess and J. Bagger Supersymmetry and Supergravity, Princeton Se- 
ries in Physics (1992) 

[12] J. Wess and B. Zumino, Phys. Lett. B 49 (1974) 52. 

[13] S. Ferrara and B. Zumino, Nucl. Phys. B 79 (1974) 413. 
A. Salam and J. Strathdee, Phys. Lett. B 51 (1974) 353. 

[14] A. Salam and J. Strathdee, Nucl. Phys. B 76 (1974) 477. 

[15] M. A. Luty and W. I. Taylor, Phys. Rev. D 53 (1996) 3399 
|arXiv:hep-th/9506098] . 



70 



M. T. Grisaru, W. Siegel and M. Rocek, Nucl. Phys. B 159 (1979) 429. 
P. C. West, Phys. Lett. B 258 (1991) 375. 

I. Jack, D. R. Jones and P. C. West, Phys. Lett. B 258 (1991) 382. 

N. Seiberg, Phys. Lett. B 318 (1993) 469 |arXiv:hep-ph/9309335| . 

M. A. Shifman and A. I. Vainshtein, Nucl. Phys. B 359 (1991) 571. 

K. A. Intriligator, R. G. Leigh and N. Seiberg, Phys. Rev. D 50 (1994) 
1092 arXiv:hep-th/9403198| . 

K. A. Intriligator, Phys. Lett. B 336 (1994) 409 |arXiv:hep-th/9407106| . 

I. Affleck, M. Dine and N. Seiberg, Nucl. Phys. B 241 (1984) 493. 
I. Affleck, M. Dine and N. Seiberg, Nucl. Phys. B 256 (1985) 557. 

G. Veneziano and S. Yankielowicz, Phys. Lett. B 113 (1982) 231. 

T. R. Taylor, G. Veneziano and S. Yankielowicz, Nucl. Phys. B 218 
(1983) 493. 

N. Dorey, T. J. Hollowood, V. V. Khoze and M. P. Mattis, Phys. Rept. 
371 (2002) 231 |arXiv:hep-th/0206063| . 

D. Finnell and P. Pouliot, Nucl. Phys. B 453 (1995) 225 
|arXiv:hep-th/9503l"T5] . 

P. Kraus and M . Shigemori, JHEP 0304 (2003) 052 
|arXiv:hep-th/0303l"04] . 

L. F. Alday and M. Cirafici, JHEP 0305 (2003) 041 
|arXiv:hep-th/0304119] . 

P. Kraus, A. V. Ryzh ov and M. Shigemori, JHEP 0305 (2003) 059 
|arXiv:hep-th/0304138l . 

M. Aganagic, K. Intriligator, C. Vafa and N. P. Warner, 
|arXiv:hep-th/030427"T| 

F. Cachazo, |arXiv:hep-th/0307063[ 

M. Matone, JHEP 0310 (2003) 068 |arXiv:hep-th/0307285| . 
K. Landsteiner and C. I. Lazaroiu, |arXiv:hep-th/0310111| 

[26] G. 't Hooft, Nucl. Phys. B 72 (1974) 461. 



71 



[27] H. Ooguri an d C. Vafa, Nucl. Phys. B 641 (2002) 3 
|arXiv:hep-th/0205297] . 



[28 
[29 

[30 
[31 
[32 
[33 

[34 
[35 
[36 

[37; 

[38 



G. Hailu and H. Georgi, arXiv:hep-th/0401101 



V. A. Novikov, M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, 
Nucl. Phys. B 229 (1983) 407. 



E. Witten, arXiv:hep-th/0302194 



P. Etingof and V. Kac, arXiv:math.qa/0305175 



E. Witten, Nucl. Phys. B 202 (1982) 253. 

A. Ceresole, G. Dall'Agata, R. D'Auria and S. Ferrara, Phys. Rev. D 
61 (2000) 066001 |arXiv:hep-th/9905226| . 

N. Seiberg, JHEP 0301 (2003) 061 |arXiv:hep-th70212225 . 



A. Gorsky, Phys. Lett. B 554 (2003) 185 |arXiv:hep-th/0210281| . 

A. Brandhuber, H. Ita, H. Nieder, Y. Oz and C. Romelsberger, 
|arXiv:hep-th /030300TJ 

R. Argurio, V. L. Campos, G. Ferretti and R. Heise, Phys. Rev. D 67 
(2003) 065005 |arXiv:hep-th/0210291| . 

M. Aganagic, A. Klemm, M. Marino and C. Vafa, 
|arXiv:hep-th/0211098| 

B. Feng, Nucl. Phys. B 661 (2003) 113 |arXiv:hep-th/0212010| . 

R. Roiban, R. Tatar and J. Walcher, Nucl. Phys. B 665 (2003) 211 
|arXiv:hep-th/0301217| . 

C. I. Lazaroiu, JHEP 0305 (2003) 044 |arXiv:hep-th7 0303008 . 

A. Klemm, K. Landsteiner, C. I. Lazaroiu and I. Runkel, JHEP 0305 
(2003) 066 |arXiv:hep-th/0303032| . 

D. Berenstein, JHEP 0306 (2003) 019 |arXiv:hep-th/0303033| . 



N. Halmagyi and V. Yasnov, arXiv:hep-th/0305134 



K. Landsteiner, C. I. Lazaroiu and R. Tatar, arXiv:hep-th/0306236 
F. Ferrari, |arXiv:hep-t h/030915l"| 



72 



[39] S. G. Naculich, H. J. Schnitzer and N. Wyllard, Nucl. Phys. B 651 
(2003) 106 |arXiv:hep-th/0211123| . 

H. Itoyama and A. Morozov, Nucl. Phys. B 657 (2003) 53 
|arXiv:hep-th/0211245] . 

S. G. Naculich, H. J. Sc hnitzer and N. Wyllard, JHEP 0301 (2003) 015 
|arXiv:hep-th/02112"54] . 

Y. Demasure and R. A. Janik, Nucl. Phys. B 661 (2003) 153 
|arXiv:hep-th/02122"l2] . 

M. Matone, Nucl. Phys. B 656 (2003) 78 |arXiv:hep-th/0212253| . 

C. h. Ahn and S. k. Nam, Phys. Rev. D 67 (2003) 105022 
|arXiv:hep-th/0301203] . 

R. Abbaspur, A. Imaanpur and S. Parvizi, JHEP 0307 (2003) 043 
|arXiv:hep-th /0302083] . 



S. G. Naculich, H. J. Schnitzer and N. Wyllard, [arXiv:hep-th/0305263 



[40] N. Dorey, T. J. Hollowood, S. Prem Kumar and A. Sinkovics, JHEP 
0211 (2002) 039 |arXiv:hep-th /0209089"] ■ 

N. Dorey, T. J. Holl owood, S. P. K umar and A. Sinkovics, JHEP 0211 
(2002) 040 |arXiv:hep-th /0209099]. 

N. Dorey, T. J. Hollow ood and S. P. Kumar, JHEP 0212 (2002) 003 
|arXiv:hep-th/0210239] . 

T. J. Hollowood, JHEP 0304 (2003) 025 |arXiv:hep-th/0212065| . 

T. Mansson, JHEP 0303 (2003) 055 | arXiv:hep-th/0302077| . 

M. Petrini, A. Tomasi ello and A. Zaffaroni, JHEP 0308 (2003) 004 
|arXiv:hep-th/030425l"] . 

T. Itoh, JHEP 0307 (2003) 005 |arXiv:hep-th/0306035| . 

[41] R. G. Leigh and M . J. Strassler, Nucl. Phys. B 447 (1995) 95 
|arXiv:hep-th/950312i"1 . 



[42] J. Polchinski and M. J. Strassler, |arXiv:hep-th/0 003136 

[43] F. Ferrari, Phys. Rev. D 67 (2003) 085013 |arXiv:hep-th/0211069| . 

F. Cachazo, N. Seibe rg and E. Witten, JHEP 0302 (2003) 042 
|arXiv:hep-th/0301006] . 

F. Ferrari, Phys. Lett. B 557 (2003) 290 |arXiv:hep-th/0301157| . 



73 



F. Cachazo, N. Seiberg and E. Witten, JHEP 0304 (2003) 018 
|arXiv:hep-th /0303207 . 

[44] T. Friedmann, Nucl. Phys. B 635 (2002) 384 |arXiv:hep-th/0203256| . 

D. Berenstein, Phys. Lett. B 552 (2003) 255 |arXiv:hep-th/0210183| . 

R. Gopakumar, JHEP 0305 (2003) 033 |arXiv:hep-th/0211100| . 

K. Ohta, JHEP 0302 (2003) 057 larXiv:hep-th/0212025| . 

C. h. Ahn and Y. Ookouchi, JHEP 0303 (2003) 010 
|arXiv:hep-th/0302150] . 

V. Balasubramanian, B. Feng, M. x. Huang and A. Naqvi, 
|arXiv:hep-th /0303065} 

T. J. Hollowood, [arXiv:hep-th/0305023 



M. Matone and L. Mazzucato, JHEP 0307 (2003) 015 
|arXiv:hep-th/03052"25] . 

M. Matone and L. Mazzucato, a rXiv:hep-th/ 0307130 

P. Svrcek, |a7x7v:hep-th/0308037l 

M. Klein and S. J. Sin, |arXiv:hep-th/0309044l 

H. Fuji and S. Mizoguchi, rarXiv:hep-th/0309049| 

C. Ahn and Y. Ookouchi, |arXiv:hep- th/0309156 , 
M. Klein and S. J. Sin, |arXiv:hep-th/0310078l 

[45] T. Eguchi and Y . Sugawara, JHEP 0305 (2003) 063 
|arXiv:hep-th/03050"50] . 

G. Bertoldi, JHEP 0306 (2003) 027 [arXiv:hep-th/0305058| . 

D. Shih, arXiv:hep-t h/030800T| 

M. A. Ganjali, |arXiv:hep-th/ 0309097 

[46] P. C. Argyres and M. R. Douglas, Nucl. Phys. B 448 (1995) 93 
|arXiv:hep-th/9505062] . 

[47] B. Feng and Y. H . He, Phys. Lett. B 562 (2003) 339 
|arXiv:hep-th/02112"34] . 

B. Feng, |arXiv:hep-th/0211202| 

B. Feng, Phys. Lett. B 572 (2003) 68 |arXiv:hep-th/0303144| . 
[48] N. Seiberg, Nucl. Phys. B 435 (1995) 129 |arXiv:hep-th /941 1 149| . 



74 



[49] A. Klemm, M. Marin o and S. Theisen, JHEP 0303 (2003) 051 
|arXiv:hep-th/0211216] . 

R. Dijkgraaf, A. Sinkovics and M. Temurhan, arXiv:hep-th/0211241 



H. Ooguri and C. Vafa, arXiv:hep-th/0302109 . 
H. Ooguri and C. Vafa, |arXiv:hep-th/0303063| 



J. R. David, E. Gava and K. S. Narain, arXiv:hep-th/0304227 



L. F. Alday, M. Cirafici, J. R. David, E. Gava and K. S. Narain, 



arXiv:hep-th/0305217 



L. F. Alday and M. Cirafici, JHEP 0309 (2003) 031 
|arXiv:hep-th /0306299] . 

F. Ardalan and N. Sadooghi, |arXiv:he p-th/0307155 . 

H. Ita, H. Nieder and Y. Qz, |arXiv:hep^tn /0309041 

R. Dijkgraaf, M. T. Grisaru, H. Ooguri, C. Vafa and D. Zanon, 
|arXiv:hep-th/031006l| 

[50] H. Fuji and Y. Ookouch i, JHEP 0212 (2002) 067 
|arXiv:hep-th/0210i48] . 

H. Ita, H. Nieder and Y. Oz, JHEP 0301 (2003) 018 
|arXiv:hep-th/02112"6T] . 

S. K. Ashok, R. Corrado, N. Halmagyi, K. D. Kennaway and C. Romels- 
berger, Phys. Rev. D 67 (2003) 086004 |arXiv:hep-th/0211291| . 

R. A. Janik and N. A. Obers, Phys. Lett. B 553 (2003) 309 
|arXiv:hep-th/0212069] . 

C. h. Ahn and S. Nam, Phys. Lett. B 562 (2003) 141 
|arXiv:hep-th/021223"T] . 

[51] Y. Ookouchi, |arXiv:hep-th/0211287| 

C. h. Ahn, Phys. Lett. B 560 (2003) 116 |arXiv:hep-th/0301011| . 

Y. Ookouchi and Y. Watabiki, Mod. Phys. Lett. A 18 (2003) 1113 
|arXiv:hep-th/03012"26] . 

C. h. Ahn, B. Feng and Y. Ookouchi, |arXiv:hep- th/0306068 , 
C. h. Ahn, B. Feng and Y. Ookouchi, |arXiv:hep- th/0307190 , 

[52] J. McGreevy, JHEP 0301 (2003) 047 |arXiv:hep-th/0211009| . 
H. Suzuki, JHEP 0303 (2003) 005 |arXiv:hep-th/0211052| . 



75 



I. Bena and R. R oiban, Phys. Lett. B 555 (2003) 117 
|arXiv:hep-th/0211075| . 

Y. Demasure and R. A. Janik, Phys. Lett. B 553 (2003) 105 
arXiv:hep-th/0211082| . 



Y. Tachikawa, arXiv:hep-th/0211189 

I. Bena, S. de Haro and R. Roiban, Nucl. Phys. B 664 (2003) 45 
|arXiv:hep-~h /0212083] . 

C. Hofman, |a7Xiv:hep-th/021 2095 , 

B. Feng, Phys. Rev. D 68 (2003) 025010 |arXiv:hep-th/0212274| . 
Y. Nakayama, JHEP 0308 (2003) 049 |arXiv:hep-th"70 306007 . 
Y. Demasure, |arXiv:hep-th/0307082| 
P. Merlatti, |arXiv:hep-th/0307115| 

B. M. Gripaios and J. F. Wheater, |arXiv:hep-th/0307176 



R. A. Janik, arXiv:hep-th/0309084. 
B. M. Gripaios, |arXiv:hep-th/0311025| 

[53] N. Seiberg, Phys. Rev. D 49 (1994) 6857 |arXiv:hep-th/9402044| . 

[54] R. Argurio, V. L. Campos, G. Ferretti and R. Heise, Phys. Lett. B 553 
(2003) 332 arXiv:hep-th/0211249| . 

I. Bena, R. Roiban and R. Tatar, |arXiv:hep-th/0211271[ 

H. Suzuki, JHEP 0303 (2003) 036 |arXiv:hep-th/0212121| . 

I. Be na, H. Murayama, R. Roiban and R. Tatar, JHEP 0305 (2003) 
049 |arXiv:hep-th^303lT5| . 



S. Corley, arXiv:hep-th/0305096 



[55] T. J. H ollowo od an d T. Kingaby, JHEP 0301 (2003) 005 
|arXiv:hep-th/0210096] . 

R. Dijkgraaf, A. Neitzke and C. Vafa, |arXiv:hep-th/02liT94"l 

S. Seki, Nucl. Phys. B 661 (2003) 257 |arXiv:hep-th/0212079| . 

S. G. Naculich, H. J. Sc hnitzer and N. Wyllard, JHEP 0308 (2003) 021 
|arXiv:hep-th /0303268 . 

R. Casero and E. Trincherini, arXiv:hep-th/0304123 . 
R. Casero and E. Trincherini, arXiv:hep-th/0307054, 



76 



[56] V. Balasubramanian, J. de Boer, B. Feng, Y. H. He, M. x. Huang, 
V. Jejjala and A. Naqvi, |arXiv:hep-t h/0212082 . 

M. Alishahiha and H. Yavartanoo, |arXiv:hep -th/0303074 , 

[57] I. Affleck, M. Dine and N. Seiberg, Phys. Lett. B 137 (1984) 187. 
I. Affleck, M. Dine and N. Seiberg, Phys. Rev. Lett. 52 (1984) 1677. 

[58] R. Argurio, G. Ferre tti and R. Heise, JHEP 0307 (2003) 044 
|arXiv:hep-th/0306l"25] . 

K. Landsteiner, C. I. Lazaroiu and R. Tatar, arXiv:hep-th/0307182 
K. Landsteiner, C. I. Lazaroiu and R. Tatar, arXiv:hep-th/0310052 



[59] Y. Tachikawa, arXiv:hep-th/0211274 



[60] J. Ambjorn and R. A. Janik, Phys. Lett. B 569 (2003) 81 
|arXiv:hep-t h /0306242] . 

[61] L. Chekhov and A. Mironov, Phys. Lett. B 552 (2003) 293 
|arXiv:hep^t h/0209085|. 

V. A. Kazakov and A. Marshakov, J. Phys. A 36 (2003) 3107 
|arXiv:hep-th/0211236| . 

H. Itoyama and A. Morozov, Phys. Lett. B 555 (2003) 287 
|arXiv:hep-th/0211259| . 

H. Itoyama and A. M orozov, Prog. Theor. Phys. 109 (2003) 433 
|arXiv:hep-th/0212032] . 

L. Chekhov, A. Marshakov, A. Mironov and D. Vasiliev, Phys. Lett. B 
562 (2003) 323 |arXiv:hep-th/0301071| . 

A. Dymarsky and V. Pestun, Phys. Rev. D 67 (2003) 125001 
|arXiv:hep-th/0301l"35l . 

H. Itoyama and A. Morozov, |arXiv:hep -th/0301136 

A. Mironov, Fortsch. Phys. 51 (2003) 781 |arXiv:hep-th7 0301196 . 

S. Aoyama and T. Masuda, arXiv:hep-th/0309232 . 



A. Alexandrov, A. Mironov and A. Morozov, arXiv:hep-th/0310113 



[62] R. Dijkgraaf, S. Gukov, V. A. Kazakov and C. Vafa, Phys. Rev. D 68 
(2003) 045007 |arXiv:hep-th/0210238|. 



77 



[63] H. Kawai, T. Kuroki and T. Morita, Nucl. Phys. B 664 (2003) 185 
|arXiv:hep-th/0303210| . 



H. Itoyama and H. Kanno, arXiv:hep-th/0304184. 

[64] R. Boels, J. de Boer, R. Duivenvoorden and J. Wijnhout, 
|arXiv:hep-th/030406T| 

M. Alishahiha and A. E. Mosaffa, JHEP 0305 (2003) 064 
|arXiv:hep-th/0304247l . 

R. Boels, J. de Boer, R. Duivenvoorden and J. Wijnhout, 
|arXiv:hep-th/0305189l 

M. Alishahiha, J. de Boer, A. E. Mosaffa and J. Wijnhout, 
|arXiv:hep-th/0308l20l 

A. Ritz, |arXiv:hep-th/0308144| 

[65] R. Dijkgraaf and C. Vafa, |arXiv:hep-th/0302011| 

T. J. Hollowood, JHEP 0303 (2003) 039 |arXiv:hep-th/0302165| . 

I. Bena and R. Roiban, |arXiv:hep-th/ 0308013 , 

T. J. Hollowood, A. Iqbal and C. Vafa, |arXiv:hep-th/0310272| 
S. Ganguli, O. J. Ganor and J. A. Gill, |arXiv:hep-th/0311042| 



78 



